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Abstract. The main stages in the evolution of a neutron star, from its birth as a 
proto-neutron star, to its old age as a cold, catalyzed configuration, are described. A 
proto-neutron star is formed in the aftermath of a successful supernova explosion and its 
evolution is dominated by neutrino diffusion. Its neutrino signal is a valuable diagnostic 
of its internal structure and composition. During its transformation from a hot, lepton- 
rich to a cold, catalyzed remnant, the possibility exists that it can collapse into a 
black hole, which abruptly terminates neutrino emissions. The essential microphysics, 
reviewed herein, that controls its evolution are the equation of state of dense matter 
and its associated neutrino opacities. Several simulations of the proto-neutron star 
evolution, involving different assumptions about the composition of dense matter, are 
described. After its evolution into a nearly isothermal neutron star a hundred or so years 
after its birth, it may be observable through its thermal emission in X-rays during its 
life in the next million years. Its surface temperature will depend upon the rapidity 
of neutrino emission processes in its core, which depends on the composition of dense 
matter and whether or not its constituents exhibit superfluidity and superconductivity. 
Observations of thermal emission offer the best hope of a determination of the radius 
of a neutron star. The implications for the underlying dense matter equation of state 
of an accurate radius determination are explored. 

1 Introduction: The Tale 

A proto-neutron star (PNS) is born in the aftermath of a successful supernova 
explosion as the stellar remnant becomes gravitationally decoupled from the 
expanding ejecta. Initially, the PNS is optically thick to neutrinos, that is, they 
are temporarily trapped within the star. The subsequent evolution of the PNS is 
dominated by ^—diffusion which first results in dclcptonization and subsequently 
in cooling. After a much longer time, photon emissions compete with neutrino 
emissions in neutron star cooling. 

In this paper, we will focus upon the essential microphysical ingredients that 
govern the macrophysical evolution of neutron stars: the equation of state (EOS) 
of dense matter and its associated neutrino opacity. Among the characteristics of 
matter that widely vary among EOS models are their relative compressibilities 
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(important in determining a neutron star's maximum mass), symmetry energies 
(important in determining the typical stellar radius and in the relative proton 
fraction) and specific heats (important in determining the local temperature). 
These characteristics play important roles in determining the matter's compo- 
sition, in particular the possible presence of additional components (such as 
hyperons, a pion or kaon condensate, or quark matter), and also significantly 
affect calculated neutrino opacities and diffusion time scales. 

The evolution of a PNS proceeds through several distinct stages ]l[D and 
with various outcomes, as shown schematically in Fig. [I]. Immediately following 
core bounce and the passage of a shock through the outer PNS's mantle, the 
star contains an unshocked, low entropy core of mass M c ~ 0.7 M in which 
neutrinos are trapped (the first schematic illustration, labelled (1) in the figure). 
The core is surrounded by a low density, high entropy (5 < s < 10) mantle 
that is both accreting matter from the outer iron core falling through the shock 
and also rapidly losing energy due to electron captures and thermal neutrino 
emission. The mantle extends up to the shock, which is temporarily stationary 
at a radius of about 200 km prior to an eventual explosion. 




cold core, worm crust isothermal (6) I0 2 <t<3x10 6 yr 

observable as X-ray 
thermal emission 

Fig. 1. The main stages of evolution of a neutron star. Shading indicates approximate 
relative temperatures. 

After a few seconds (stage 2), accretion becomes less important if the super- 
nova is successful and the shock lifts off the stellar envelope. Extensive neutrino 
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losses and deleptonization will have led to a loss of lepton pressure and the col- 
lapse of the mantle. If enough accretion occurs, however, the star's mass could 
increase beyond the maximum mass capable of being supported by the hot, 
lepton-rich matter. If this occurs, the remnant collapses to form a black hole 
and its neutrino emission is believed to quickly cease 

Neutrino diffusion deleptonizes the core on time scales of 10-15 s (stage 
3). The diffusion of high-energy (200-300 MeV) neutrinos from the core to the 
surface where they escape as low-energy neutrinos (10-20 MeV) generates a 
large amount of heat within the star (a process akin to joule heating). The 
core's entropy approximately doubles, producing temperatures in the range of 
30-60 MeV, during this time, even as neutrinos continue to be prodiguously 
emitted from the stars effective surface, known as the ^—sphere. 

Strange matter, in the form of hyperons, a Bose condensate, or quark mat- 
ter, which is suppressed to extremely large densities when neutrinos are trapped 
in matter, could appear at the end of the deleptonization. The appearance of 
strange matter leads to a decrease in the theoretical maximum mass that mat- 
ter is capable of supporting, leading to another possibility for black hole for- 
mation B. This would occur if the PNS's mass, which must be less than the 
maximum mass of hot, lepton-rich matter (or else a black hole would already 
have formed), is greater than the maximum mass of hot, lepton-poor matter. 
However, if strangeness does not appear, the theoretical maximum mass instead 
increases during deleptonization and the appearance of a black hole would be 
unlikely unless accretion in this stage remains significant. 

The PNS is now lepton-poor, but it is still hot. While the star has zero net 
neutrino number, thermally produced neutrino pairs of all flavors are abundant 
and dominate the emission. Neutrino diffusion continues to cool the star, but the 
average neutrino energy decreases, and the neutrino mean free path increases. 
After approximately 50 seconds (stage 4), the mean free path becomes compa- 
rable to the stellar radius, and star finally becomes transparent to neutrinos. 
Since the threshold density for the appearance of strange matter decreases with 
decreasing temperature, a delayed collapse to a black hole is still possible during 
this epoch. 

Neutrino observations from a galactic supernova will illuminate these stages. 
The observables will constrain time scales for deleptonization and cooling and 
the star's binding energy. Dimcnsionally, diffusion time scales are proportional to 
i? 2 (cA) _1 , where R is the star's radius and A is the effective neutrino mean free 
path. This generic relation illustrates how both the EOS and the composition, 
which determine both R and A, influence evolutionary time scales. The total 
binding energy, which is primarily a function of stellar mass and radius (Lattimer 
& Prakash Q), should be one of the most accurately measured quantities from 
neutrino observatories. Currently, Super-Kamiokande and SNO are capable of 
detecting thousands of neutrinos from a galactic supernova (distance less than 
10 kpc). Exciting possibilities lie ahead with many other existing and planned 
new facilities Q. 
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Following the onset of neutrino transparency, the core continues to cool by 
neutrino emission, but the star's crust remains warm and cools less quickly. The 
crust is an insulating blanket which prevents the star from coming to complete 
thermal equilibrium and keeps the surface relatively warm (Tk3x 10 6 K) for 
up to 100 years (stage 5). This timescale is primarily sensitive to the neutron 
star's radius and the thermal conductivity of the mantle H, as can be noted from 
the approximate diffusive relationship r oc AR 2 /X, where AR is the thickness 
of the crust. If the rapid decrease in the star's surface temperature predicted to 
occur when thermal equilibrium is ultimately achieved (see Fig. [l6] in Section 
6), a valuable constraint on the thickness of the crust, and hence the neutron 
star radius, could be obtained. The temperature of the surface after the interior 
of the star becomes isothermal (stage 6) is determined by the rate of neutrino 
emission in the star's core. The magnitude of the rate is primarily determined by 
the question of whether or not one or more of the so-called direct Urea processes 
can occur. The basic Urea process 

n — > p + e~ + D e ; p—tn + e + + v e (1) 

operates even in degenerate matter because at finite temperature some of the 
nucleons are in excited states. In addition, direct Urea process involving hyper- 
ons, Bose condensates and quarks are also possible. In general, the direct Urea 
rate is proportional to T 4 , and is so large that the surface temperatures fall to 
just a few times 10 5 K, which becomes very difficult to observe in X-rays ex- 
cept for very nearby stars. A relatively high surface temperature, closer to 10 6 
K, will persist, however, if an Urea process can only occur indirectly with the 
participation of a spectator nucleon - the modified Urea process, which in the 
case of nucleons is 

n + (n,p) — ► p + (n,p) + e~ + i> e ; p + (n,p) — > n + {n 7 p) + e + + v e , (2) 

and leads to the so-called standard cooling scenario. 

However, there are two circumstances that could prevent the direct Urea 
process from occuring. First, if the composition of the matter is such that the 
momentum triangle involving the non-neutrino particles cannot be closed, mo- 
mentum conservation disallows this process. This occurs, in the case of n,p,e, 
for example, if the p and e~ abundances, which must be equal, are less than 
1/8 the n abundance. This would be the case if the nuclear symmetry energy 
has a relatively weak density dependence. In addition, direct Urea processes in- 
volving hyperons, a Bose condensate, or quarks would not occur, of course, if 
they are not present. Second, direct Urea processes are suppressed if one of the 
reactants becomes superfluid. In this case, when the core temperature falls be- 
low the superfluid's critical temperature, the rapid cooling is terminated. In the 
case of a superfluid, the core cooling, and therefore the surface temperature, will 
be intermediate between those predicted by standard and rapid cooling models. 
Neutrino emission continues to dominate until neutron stars are approximately 1 
million years old, at which point photon cooling from the surface begins to dom- 
inate. Unless the interiors cool very rapidly, X-ray emissions from stars remain 
relatively high until the photon cooling epoch. 
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Several neutron stars have been suggested to have observable thermal emis- 
sions in X-rays. In addition, the nearby neutron star RX J185635-3754, which 
is the closest known neutron star with a distance of approximately 60 pc, has 
detectable UV and optical thermal emissions as well. Such objects represent the 
best chance for measuring a neutron star's radius, especially if the redshift of 
spectral lines can be determined. Just-launched or proposed X-ray satellites, 
such as Chandra and XMM offer abundant prospects of observations of photon 
observations of neutron stars. 

The organization of this article is as follows. Section 2 contains a summary of 
the basic equations of evolution for proto-neutron stars, including a discussion of 
the equilibrium diffusion approximation. Section 3 details the equation of state 
of dense matter, taking into account the possibility that neutrinos are trapped 
in the matter. The possibilities of hyperon-rich matter, kaon condensation and 
strange quark matter are also discussed here. In the event of a transition to 
matter containing kaons or quarks, we also consider the possibility that matter 
could be inhomogeneous with droplets of the strange matter embedded in normal 
matter. Neutrino-matter interactions are considered in Section 4, which includes 
discussions of the effects of composition, in-mcdium dispersion relations, and 
correlations, in both homogeneous and inhomogeneous phases. In Section 5, we 
present several simulations of the evolution of proto-neutron stars with different 
assumptions about the composition of dense matter and highlighting the role of 
the neutrino opacities. Focus is placed upon predicted neutrino signals and the 
differences anticipated for varying assumptions about the matter's composition. 
Section 6 describes the long-term cooling epoch, with a special emphasis on the 
role of direct Urea processes and superfluidity and a comparison with observa- 
tions. A discussion of the possibility of detecting superfluidity, including quark 
color superfluidity, is included. In Section 7, the dependence of the structure of 
neutron stars on the underlying dense matter equation of state is explored. The 
relation between the matter's symmetry energy and the radii of neutron stars 
is highlighted. In addition, the moments of inertia and binding energies of neu- 
tron stars are discussed, and observational constraints on the mass and radius 
of the Vela pulsar from these considerations are elaborated. Section 8 contains 
our outlook. 



2 Short-Term Neutrino Cooling: The First Minutes 

The cooling of PNSs can be divided into two main regimes: the short-term, 
lasting perhaps one minute, during which the potential to observe the neutrino 
signal in terrestrial detectors exists, and the longer term period, lasting perhaps 
one million years, in which neutrino emissions dominate the cooling but the 
star is observable only through its thermal, photonic, emissions. This section 
summarizes the evolution equations relevant for the short-term Kelvin- Hclmholtz 
phase and the estimation of its neutrino signature. 
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2.1 PNS Evolution Equations 

The equations that govern the transport of energy and lepton number are ob- 
tained from the Boltzmann equation for massless particles [|ijP,|9],|l0t . We will 
focus on the non-magnetic, spherically symmetric situation. For the PNS prob- 
lem, fluid velocities are small enough so that hydrostatic equilibrium is nearly 
fulfilled. Under these conditions, the neutrino transport equations in a stationary 
metric 

ds 2 = -e^dt 2 + e 2A dr 2 + r 2 d9 2 + r 2 sin 2 9 d<P 2 (3) 
d{N v /n B ) , d{eH<Kr 2 F u ) _ ,S N 



arc: 



dt da n B 



(4) 



W , p d(l/n B ) _^ d{e 2 Hnr 2 H v ) _ ^S E 

v r e — e , {0 

at at da n B 

where n B is the baryon number density and a is the enclosed baryon number 
inside a sphere of radius r. The quantities N v , P„, and Sn are the number 
density, number flux and number source term, respectively, while J„, H v , P„, 
and Se are the neutrino energy density, energy flux, pressure, and the energy 
source term, respectively. 

In the absence of accretion, the enclosed baryon number a is a convenient 
Lagrangian variable. The equations to be solved split naturally into a trans- 
port part, which has a strong time dependence, and a structure part, in which 
evolution is much slower. Explicitly, the structure equations are 

dr 1 dm p 

da Airr 2 n B e A ' da n B e A 

^ = ^—( m + 4nr s P) , ^ = -(„ + P) T ^—(m + 4*r*P) . (7) 
da iirr^ns da Aitr^ns 

The quantities m (enclosed gravitational mass), p (mass-energy density), and 
P (pressure) include contributions from the leptons. To obtain the equations 
employed in the transport, Eq. (^) may be combined with the corresponding 
equation for the electron fraction 

^=-e^ (8) 
dt n B w 



to obtain 

dY L ,d(e^^r 2 F v ) 



0. (9) 



dt da 
Similarly, Eq. (§) may be combined with the matter energy equation 

dt dt n B 
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where U is the specific internal energy and use of the first law of thermodynamics 
yields 

6 at + ^ dt + da ~ u ' 1 ' 

where s is the entropy per baryon. 



2.2 The Equilibrium Diffusion Approximation 



At high density and for temperatures above several MeV, the source terms in the 
Boltzmann equation are sufficiently strong to ensure that neutrinos are in ther- 
mal and chemical equilibrium with the ambient matter. Thus, the neutrino dis- 
tribution function in these regions is both nearly Fermi-Dirac and isotropic. We 
can approximate the distribution function as an expansion in terms of Legendre 
polynomials to O(fi), which is known as the diffusion approximation. Explicitly, 



/(w,M)=/o(w)+f*/i(w), fo = [l + ei^)]- 1 



(12) 



where fo is the Fermi-Dirac distribution function at equilibrium (T = T mat , 
[i v — fif, q ), with u> and [i v being the neutrino energy and chemical potential, 
respectively. The main goal is to obtain a relation for f\ in terms of fo- In the 
diffusion approximation, one obtains |l(| 



h 



-D{ui) 



ydfo 



{dcj> df 



dr dr du> 

The explicit form of the diffusion coefficient D appearing above is given by 



D{w) 



J + 



(13) 



(14) 



The quantity j — j a + j s , where j a is the emissivity and j s is the scattering 
contribution to the source term. The absorptivity is denoted by A a and fcf is the 
scattering contribution to the transport opacity. Substituting 



dfo 
dr 



T 



dr 



oj dT 
T~dr~ 



dfo 
duj 



where r\ v 



H v /T is the neutrino degeneracy parameter, in Eq. 

u d(Te^) 



fi 



-D(uj)e 



dr Te<t> dr 



dfo 
duj 



(15) 

3|), we obtain 
(16) 



Thus, the energy-integrated lepton and energy fluxes are 



F v = — 
H v = - 



e -V 



6ir 2 



-Cj>rp3 



6tt 2 



D. 



diTe^) 

dr 
d(Te*) 

dr 



(Te^ 2 g 



{Te*)D- 



di] 
dr 



(17) 
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The coefficients D 2 , D 3 , and D4 are related to the energy-dependent diffusion 
coefficient D{uj) through 



where x = ui/T. These diffusion coefficients depend only on the microphysics of 
the neutrino-matter interactions (see §4 for details). The fluxes appearing in the 
above equations are for one particle species. To include all six neutrino types, 
we redefine the diffusion coefficients in Eq. (|17|) : 



D 2 = + Dl" , D 3 = Dp ~ , D 4 = + D»° + AD V A » . (19) 



2.3 Neutrino Luminosities 

A fair representation of the signal in a terrestrial detector can be found from 
the time dependence of the total neutrino luminosity and average neutrino en- 
ergy together with an assumption of a Fermi-Dirac spectrum with zero chemical 
potential. We will return to discuss the improvements necessary to obtain more 
accurate information about the spectra. 

The total neutrino luminosity is the time rate of change of the star's gravita- 
tional mass, and is therefore primarily a global property of the evolution. This 
luminosity, due to energy conservation, must also equal 



at the edge of the star. This relation serves as a test of energy conservation, at 
least for all times greater than about 5 ms, when the star comes into radiative 
equilibrium. For times greater than about 5 ms, initial transients become quite 
small and the predicted luminosities should be relatively accurate compared to 
full transport simulation. Estimate of the average energy of neutrinos is made 
from the temperature T v of the matter at the neutrinosphere R u , defined to 
be the location in the star where the flux factor = 0.25. However, since 
the spectrum may not be Fermi-Dirac at the neutrinosphere, a diffusion scheme 
cannot give a very precise value for the average energy. We use the average energy 
< E u 3T„, where T v is a mass average in the outermost zone. Because it is 
a globally determined quantity, the luminosity L v is necessarily more accurately 
determined than either R v or T v . 

3 The Equation of State of Neutrino Trapped Matter 

The rationale for considering different possibilities for the composition of dense 
matter is largely due to the fact that QCD at finite baryon density remains 
unsolved. Effective QCD based models have raised intriguing possibilites con- 
cerning the composition of dense matter including the presence of hyperons, 
pion or kaon condensates, and quark matter (see Q for extensive references). It 




(18) 



(20) 
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is also important to have predictions for the plain- vanilla case of nucleons alone. 
The contrast can be dramatic, since additional components offer the possibility 
of BH formation during the evolution of a PNS. In what follows, the symbols np 
refer to matter with nucleons alone, npH to matter including hyperons, npK to 
matter with nucleons and kaons, and npQ to matter with nucleons and quarks. 
In all cases, leptons in beta equilibrium are included. 

3.1 Matter with Nucleons and Hyperons 

The masses and radii of neutron stars depend upon the matters' compressibility, 
the composition of matter at high density, and the nuclear symmetry energy 
(e.g., [Q). In the PNS problem, the finite temperature aspects of the EOS also 
play an important role. During the early evolution the entropy in the central 
regions is moderately high, s ~ 1 — 2 (in units of Boltzmann's constant), which 
correspond to temperatures in the range T = 20 — 50 MeV. These features may 
be explored by employing a finite temperature field-theoretical model in which 
the interactions between baryons are mediated by the exchange of a, ui, and p 
mesons^. The hadronic Lagrangian density is given by |Q 

L H = ^2 B^-ij^dfj - gwi^uJ^ - s^b^ • t - M, + g„iu)Bi 

i 

- \w, v W^ + \ml^ - \b, v B^ + \m%b» 

+ l -d^a~ l -mla 2 -U{o) (21) 

Here, B are the Dirac spinors for baryons and t is the isospin operator. The sums 
include baryons i = n,p, A, S, and S. The field strength tensors for the tu and 
p mesons are W M „ = d^ v — <9„w M and B M „ = <9 M b„ — d^h^, respectively. The 
potential U(a) represents the self-interactions of the scalar field and is taken to 
be of the form Jil| 

U{a) = hjM n (g aN af + ^c(g aN a) 4 . (22) 

The partition function Zh for the hadronic degrees of freedom is evaluated in the 
mean field approximation. The total partition function Z to tai = ZhZl, where Zl 
is the standard noninteracting partition function of the leptons. Using Z to tai-, the 
thermodynamic quantities can be obtained in the standard way. The additional 
conditions needed to obtain a solution are provided by the charge neutrality 

1 Note that the couplings in these models may be chosen to reproduce the results 
of numerically more intensive microscopic potential models, such as that of Akmal 
and Pandharipande , so that the gross features of the zero temperature thermo- 
dynamics can be reproduced. Additional advantages to this approach are that the 
effects of finite temperature and arbitrary proton fraction may be incorporated more 
easily. 
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requirement, and, when neutrinos are trapped, the set of equilibrium chemical 
potential relations required by the general condition 

Mi = &iMn-ft(w ( 23 ) 

where bi is the baryon number of particle i and qi is its charge. The introduc- 
tion of additional variables, the neutrino chemical potentials, requires additional 
constraints, which we supply by fixing the lepton fractions, Yn, appropriate for 
conditions prevailing in the evolution of the PNS. In addition to models con- 
taining only nucleonic degrees of freedom (GMlnp & GM3np) we investigate 
models that allow for the presence of hyperons (GMlnpH & GM3npHV. For the 
determination of the various coupling constants appearing in Zh see |4j . 

The lepton chemical potentials influence the deleptonization epoch. For np 
models a lower nuclear symmetry energy favors a larger v e fraction and has little 
effect on the e~ fraction at Yl — 0.4. Models with hyperons lead to significantly 
larger and lower \x e , both of which influence the diffusion of electron neutri- 
nos. The electron chemical potentials in neutrino free matter are reduced to a 
greater extent by changes in composition and symmetry energy as there are no 
neutrinos to compensate for changes in p, = \L n — [i v . 



3.2 Matter with a Kaon Condensate 

The contents of this section are extracted from Pons et al. Q|. For the kaon 
sector, we use a Lagrangian that contains the usual kinetic energy and mass 
terms along with the meson interactions Kaons are coupled to the meson 
fields through minimal coupling; specifically, 

L K = V*K + V^K~ - m^K+K- , (24) 

where the vector fields are coupled via the standard form 

T>n = + ig^K^p. + igpKl^bf, ■ t (25) 

and m* K = rriK — ^oK® is the effective kaon mass. 

In the mean field approach, the thermodynamic potential per unit volume in 
the kaon sector is |l4| 

^ = !W«-(m + *o) 2 ) 



T 



d 3 p 



ln(l-e-^ w -^)+\n(l-e-^ ++ ^)] , (26) 



(2tt) 3 
o 

where Xq = g^K^o + gpxbo, the Bose occupation probability 



f B (x) = (e 13 * - I)" 1 , w ± = yV 2 + m^ 2 ± X a , f = 93 MeV is the pion decay 
constant and the condensate amplitude, 9, can be found by extremization of 
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the partition function. This yields the solution 8 = (no condensate) or, if a 
condensate exists, the equation 

m *i< = + X (27) 

where \ik is the kaon chemical potential. In beta-stable stellar matter the con- 
ditions of charge neutrality 

^ qBns —n e — n K = (28) 

B 

and chemical equilibrium 

Hi = biH n - qi(m - n„ t ) (29) 
= A*n - Hp (30) 

are also fulfilled. 

The kaon condensate is assumed to appear by forming a mixed phase with 
the baryons satisfying Gibbs' rules for phase equilibrium Jig] . Matter in this 
mixed phase is in mechanical, thermal and chemical equilibrium, so that 

p '=p", T' = T", n{ = nl 1 , (31) 

where the superscripts I and II denote the nucleon and kaon condensate phases, 
respectively. The conditions of global charge neutrality and baryon number con- 
servation are imposed through the relations 

Xq 1 + (1 - x)q H =0 

X n B + (1 - xWb 1 = n B , (32) 

where x denotes the volume fraction of nucleonic phase, q the charge density, 
and n b the baryon density. We ignore the fact that the phase with the smallest 
volume fraction forms finite-size droplets; in general, this would tend to decrease 
the extent of the mixed phase region. Further general consequences of imposing 
Gibbs' rules in a multicomponent system are that the pressure varies continu- 
ously with density in the mixed phase and that the charge densities must have 
opposite signs in the two phases to satisfy global charge neutrality. We note, 
however, that not all choices of nucleon-nucleon and kaon-nucleon interactions 
permit the Gibbs' rules to be satisfied (for an example of such an exception, see 
| fl4| ) . The models chosen in this work do allow the Gibbs' rules to be fulfilled at 
zero and finite temperatures and in the presence of trapped neutrinos. 

The nucleon-meson couplings are determined by adjusting them to reproduce 
the properties of equilibrium nucleonic matter at T = 0. We use the numerical 
values used by [ fl3"|| , i.e., equilibrium density no = 0.153 fm -3 , equilibrium energy 
per particle of symmetric nuclear E/A = —16.3 MeV, effective mass M* = 
0.78M, compression modulus Kq = 240 MeV, and symmetry energy a sym = 32.5 
MeV. These values yield the coupling constants g a /m a = 3.1507 fm, g^/m^ = 
2.1954 fm, g p /m p = 2.1888, b = 0.008659, and c = -0.002421. 
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The kaon-meson couplings g a x and g^K are related to the magnitude of 
the kaon optical potential Uk at the saturation density no of isospin symmetric 
nuclear matter: 

U K (n ) = -gaK<y{na) - g u K^o(no). (33) 

Fits to kaonic atom data have yielded values in the range —(50 — 200) MeV 
[ ^7l|l^ , [l9| , ^0|j2l|l . We use g U K — 3wJv/3 and g P K — g P N /2 on the basis of simple 
quark and isospin counting. Given the uncertainty in the magnitude of \Uk\, 
consequences for several values of \Uk\ were explored in [[u]. Moderate values 
of \Uk\ generally produce a second order phase transition and, therefore, lead 
to moderate effects on the gross properties of stellar structure. Values in excess 
of 100 MeV were found necessary for a first order phase transition to occur; in 
this case kaon condensation occurs at a relatively low density with an extended 
mixed phase region, which leads to more pronounced effects on the structure due 
to a significant softening of the EOS. 

1.2 
1.0 

0.8 
J 0.6 
0.4 
0.2 
0.0 

0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 

Fig. 2. Phase boundaries between pure nucleonic matter (I), pure kaon condensed 
matter (II) and a mixed phase (III) in the Yl-ub plane for Uk = —100 MeV (left 
panel) and Uk ~ —120 MeV (right panel). The solid line corresponds to s = and the 
dashed line to s = 1. The dashed-dotted line shows the baryon density as a function of 
the lepton fraction for s — 0, neutrino-free (Yl = Y e ) matter 

The phase boundaries of the different phases are displayed in Fig. || in a 
Y L -n B plane for an optical potential U K of -100 MeV (left) and -120 MeV 
(right), respectively. The nucleonic phase, the pure kaon matter phase, and the 
mixed phase are labelled I, II, and III, respectively. Solid lines mark the phase 
transition at zero temperature and dashed lines mark the phase transition at an 
entropy per baryon of s = 1. Note that finite entropy effects are small and do 
not affect significantly the phase transition density. The dash-dotted line shows 
the electron fraction Y e as a function of density in cold, catalyzed matter (for 
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which Yl = Ye), which is the final evolutionary state. The region to the left 
of this line corresponds to negative neutrino chemical potentials and cannot be 
reached during normal evolutions. The solid and dashed lines, which separate 
the pure phases from the mixed phase, vary roughly linearly with the lepton 
fraction. Also notice the large, and nearly constant, densities of the boundary 
between the mixed phase III and the pure kaon phase II. These densities, for 
the cases shown, lie above the central densities of the maximum mass stars, so 
that region II does not generally exist in proto-neutron stars (see (l4|). The 
effect of increasing the lepton number is to reduce the size of the mixed phase 
(which in fact shrinks to become a second order phase transition for Yl > 0.4 
and Uk = —100 MeV) and to shift the critical density to higher densities. A 
similar effect is produced by decreasing the magnitude of the optical potential. 

3.3 Matter with Quarks 

The discussion in this section follows the work of Steiner, Prakash & Lattimer 
[ p2| . The thermodynamic potential of the quark phase is D = J?fg + ^int, where 




denotes the Fermi gas contribution arising from quarks. We consider three fla- 
vors, i = u,d,s and three colors, N c = 3 of quarks. The distribution functions 
of fermions and anti-fermions are /j = [1 + exp(/3(Ei — and 
fi = [1 + exp(/3(-E7, + ju,i))] , where Ei and fa are the single particle energy and 
chemical potential, respectively, of quark species i. To explore the sensitivity 
of the quark model, we contrast the results of the MIT bag and the Nambu 
Jona-Lasinio (henceforth NJL) models for J?i n t- 

In the MIT bag model, the Fermi gas contribution is calculated using current, 
as opposed to dynamical, quark masses. We will restrict ourselves to the simplest 
bag model and keep only the constant cavity pressure term. The results are 
qualitatively similar to what is obtained by including perturbative corrections, 
if the bag constant B is slightly altered (2^] . 

Several features of the Lagrangian of Quantum Chromo-Dynamics (QCD), 
including the spontaneous breakdown of chiral symmetry, are exhibited by the 
Nambu Jona-Lasinio (NJL) model, which shares many symmetries with QCD. 
In its commonly used form, the NJL Lagrangian reads 

8 

C = q(ip - mo)q + G*y^[ (qA fc <?) 2 + (gijsXkq) 2 } 

fc=0 

- K[detf(q(l + 75 )q)+det f (q(l- 75 )q)} . (35) 

The determinant operates over flavor space, rh is the 3x3 diagonal current 
quark mass matrix, A^ represents the 8 generators of SU(3), and Ao is pro- 
portional to the identity matrix. The four-fermion interactions stem from the 
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original formulation of this model |24[ , while the flavor mixing, determincntal 
interaction is added to break {7a (1) symmetry [ p5[ . Since the coupling constants 
G and K are dimensionful, the quantum theory is non-renormalizable. Therefore, 
an ultraviolet cutoff A is imposed, and results are considered meaningful only 
if the quark Fermi momenta are well below this cutoff. The coupling constants 
G and K, the strange quark mass m S fi, and the three-momentum ultraviolet 
cutoff parameter A, are fixed by fitting the experimental values of /„•, m w , mx 
and m v >. We use the values of Ref. |§, namely A = 602.3 MeV, GA 2 = 1.835, 
KA 5 — 12.36, and m ^ s — 140.7 MeV, obtained using mo,u = mo,d = 5.5 MeV. 
The subscript "0" denotes current quark masses. Results of the gross properties 
of PNSs obained by the alternative parameter sets of Refs. (2^] and ^8| are 
similar to the results quoted below. 

In the mean field approximation at finite temperature and at finite baryon 
density, the thermodynamic potential due to interactions among quarks is given 

^Int oAr f d3 P 



2G{qiq t ) 2 - AK(q u q u ) (qaqd) (q s q s ) ■ (36) 



In both Eqs. (34) and (|36|) for the NJL model, the quark masses are dynami- 
cally generated as solutions of the gap equation obtained by requiring that the 
potential be stationary with respect to variations in the quark condensate (qiQi): 

m l = m 0<i - 4G(q~iqi) + 2K {qjqj)(q k q k ) , (37) 

(qi, qj, gx.) representing any permutation of (u, d, s). The quark condensate 
and the quark number density nj = (qjqi) are given by: 



{qiOi) =-2N c I £^3^ [1 -h-h 



d 3 p rrii 



n t = (q}q t ) = 2N C [ [f t - /,] . (38) 

J (27T) 

A comparison between the MIT bag and NJL models is facilitated by defining 



an effective bag pressure in the NJL model to be |29 B e // = f^int/V — B 
with BqV = fii n t\n u =n d =n 3 =Q & constant value which makes the vacuum energy 
density zero. In this way, the thermodynamic potential can be expressed as 
Q = B e f fV+SlpQ which is to be compared to the MIT bag result Q = BV '+J7fg- 
Note, however, that Qfg m the NJL model is calculated using the dynamical 
quark masses from Eq. (|37]) . 

The temperature as a function of baryon density for fixed entropy and net lep- 
ton concentration is presented in Fig. |^, which compares the cases (s = 1, Yl s = 
0.4) and (s = 2, Y Vi! = 0) both including and ignoring quarks. The temperature 
for a multicomponent system in a pure phase can be analyzed with the relation 
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Fig. 3. Temperature versus density in units of no for two PNS evolutionary snapshots. 
The upper (lower) panel displays results for the NJL (MIT bag) Lagrangian. The 
parameters £ = £ = in the Miiller-Serot (MS) hadronic Lagrangian [B0| are chosen. 
Results are compared for matter containing only nucleons (np), nucleons plus hyperons 
(npH), nucleons plus quarks (npQ) and nucleons, hyperons and quarks (npHQ). Bold 
curves indicate the mixed phase region 



for degenerate Fermi particles 




where m* and are the effective mass and the Fermi momentum of component 
i, respectively. This formula is quite accurate since the hadron and quark Fermi 
energies are large compared to the temperature. The introduction of hyperons 
or quarks lowers the Fermi energies of the nucleons and simultaneously increases 
the specific heat of the matter, simply because there are more components. In 
the case of quarks, a further increase, which is just as significant, occurs due to 
the fact that quarks are rather more relativistic than hadrons. The combined 
effects for quarks are so large that, in the case Mq = 0.6M shown in Fig. ||, an 
actual reduction of temperature with increasing density occurs along an adiabat. 
The effect is not necessarily as dramatic for other choices of Mq , but nevertheless 
indicates that the temperature will be smaller in a PNS containing quarks than 
in stars without quarks. The large reduction in temperature might also influence 
neutrino opacities, which are generally proportional to T 2 . However, the presence 
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of droplet-like structures in the mixed phase, not considered here, will modify 
the specific heat. In addition, these structures may dominate the opacity in 
the mixed phase However, a PNS simulation is necessary to consistently 
evaluate the thermal evolution, since the smaller pressure of quark-containing 
matter would tend to increase the star's density and would oppose this effect. 



1.2 
I 



- s = 

1 s =l MS-NJL J 

I s = a npQ : 

II F 




7 MS MIT \ n - 
I B = 200 MeV/frrr 3 ; 
7 npQ 7 




- 7 




I 

1 — 1 1 1 III 1 — Mllll" 


7 I 1 

: — 1 — 1 1 1 1 " — i — 1 1 1 1 1 r 


1 MS-NJL u 
: npHQ 


— ^— — ~ 


MS -MIT \l ~ 
npHQ 






n c 




I 

r 


I - 

r 



10- 



io- 

Y,,. 



Fig. 4. The phase diagram of the quark-hadron transition in the baryon number density 
- neutrino concentration plane for three representative snapshots during the evolution 
of a proto-neutron star. The left (right) panels are for the NJL (MIT bag) quark 
EOS, and hyperons are (are not) included in the bottom (top) panels. The parameters 
£ = £ = in the Miiller-Serot (MS) hadronic Lagrangian are chosen. The lower- and 
upper-density boundaries of the mixed phase are indicated by bold curves. The central 
densities of maximum mass configurations are shown by thin curves 



This last point is highlighted in Fig. |4| which shows phase diagrams for the 
mixed phase in the baryon density-neutrino fraction plane. The upper and lower 
boundaries of the mixed phase region are displayed as bold lines, while the 
central densities of the maximum mass configurations are shown as light lines. 
In no case, for either quark model and whether or not hyperons are included, 
are pure quark stars possible. The high-density phase boundaries are always well 
above the central densities. While in the optimum case, in which the parameters 
of both the hadronic and quark EOSs are fine-tuned, it is possible for a pure 
quark core to form if B < 150 MeV fm -3 , the maximum mass decreases below 
1.44 Mq if B < 145 MeV fm~ 3 . This narrow window, which further decreases or 
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disappears completely if the hadronic EOS is altered, suggests that pure quark 
configurations may be unlikely. 

3.4 Inhomogeneous Phases 

It is widely believed that some type of phase transition will occur in nuclear 
matter at high densities. For example, a transition to deconfined quark matter 
should exist at sufficiently high density, and at lower densities, a first-order 
transition to a Bose condensate phase might exist. Such phase transitions are 
expected to soften the equation of state, leading to changes in the mass-radius 
relation and lowering the maximum mass. Phase transitions can also influence 
transport and weak interaction rates in matter. 

Glendcnning has shown that, due to the existence of two conserved charges 
(baryon number and charge) instead of just one, first order phase transitions can 
lead to a large mixed phase region in the neutron star interior (32). The mixed 
phase consists of high baryon density, negatively charged, matter coexisting with 
lower density, positively charged, baryonic matter. The situation is entirely anal- 
ogous to the well-known situation involving the mixed phase consisting of nuclei 
and a surrounding nucleonic vapor that occurs below nuclear saturation density 
[ p3| . The occurence of a mixed phase, as opposed to a Maxwell construction, 
results in a wider transition in which bulk thermodynamic properties such as 
pressure vary less rapidly but are softer over a wider density range. In addi- 
tion, the propagation of neutrinos whose wavelength is greater than the typical 
droplet size and less than the inter-droplet spacing will be greatly affected by the 
heterogeneity of the mixed phase, as a consequence of the coherent scattering of 
neutrinos from the matter in the droplet. The thermodynamics and the effect on 
neutron star structure of two situations have been studied in some detail: first 
order kaon condensation p^| , p4p5| ] and the quark- hadron transition |3(||37j . 

4 Neutrino-Matter Interaction Rates 

One of the important microphysical inputs in PNS simulations is the neutrino 
opacity at supra-nuclear density |l|, ^ , ^9po|JIl| . |42| Although it was realized over 
a decade ago that the effects due to degeneracy and strong interactions sig- 
nificantly alter the neutrino mean free paths, it is only recently that detailed 
calculations have become available [^l]^|,||,||,|6|j4^,||,|9| . The scattering and 
absorption reactions that contribute to the neutrino opacity are 

v e + B —> e~ + B' , v e + B-> e+ + B' , (40) 
v x + B — ► v x + B' , v x + — > v x + e~ , (41) 

where the scattering reactions are common to all neutrino species and the dom- 
inant source of opacity for the electron neutrinos is due to the charged reaction. 
The weak interaction rates in hot and dense matter are modified due to many 
in-medium effects. The most important of these are: 



18 Madappa Prakash et al. 



(1) Composition: The neutrino mean free paths depend sensitively on the compo- 
sition which is sensitive to the nature of strong interactions. First, the different 
degeneracies of the different Fermions determines the single-pair response due to 
Pauli blocking. Second, neutrinos couple differently to different baryonic species; 
consequently, the net rates will depend on the individual concentrations. 

(2) In-medium dispersion relations: At high density, the single-particle spectra 
are significantly modified from their noninteracting forms due to effects of strong 
interactions. Interacting matter features smaller effective baryon masses and en- 
ergy shifts relative to non-interacting matter. 

(3) Correlations: Repulsive particle-hole interactions and Coulomb interactions 
generally result in a screened dielectric response and also lead to collective ex- 
citations in matter. These effects may be calculated using the Random Phase 
Approximation (RPA) , in which ring diagrams are summed to all orders. Model 
calculations Jl3[i5 4?|. ^o| , p r lp^ , |53] ,p4 55 indicate that at high density the neu- 



trino cross sections are suppressed relative to the case in which these effects 
are ignored. In addition, these correlations enhance the average energy transfer 
in neutrino-nucleon collisions. Improvements in determining the many-body dy- 
namic form factor and assessing the role of particle-particle interactions in dense 
matter at finite temperature are necessary before the full effects of many-body 
correlations may be ascertained. 

The relative importance of the various effects described above on neutrino 
transport is only beginning to be studied systematically. As a first step, we 
will focus on effects due to modifications (1) through (3) above. To see how 
this is accomplished, we start with a general expression for the differential cross 



section |56 57 



1 d 3 a G% E 3 



V d 2 Q 3 dE 3 128tt 2 Ei 



, qo- (ms - im) 
1 - ex P 1 f 



x (1 - f 3 (E 3 )) Im (L^7T« ) , (42) 

where the incoming neutrino energy is E\ and the outgoing electron energy is 
E 3 . The factor [1 — exp((— qo — /i2 + fi^/T)]^ 1 maintains detailed balance, for 
particles labeled '2' and '4' which are in thermal equilibrium at temperature T 
and in chemical equilibrium with chemical potentials fi2 and respectively. 
The final state blocking of the outgoing lepton is accounted for by the Pauli 
blocking factor (1 — f3(E 3 )). The lepton tensor L a p is given by 

L a > 3 = 8[2fc a fc /3 + (k ■ q)g af} - (fcV + q a k fj ) =p it af3 ^k^q u ] (43) 

The target particle retarded polarization tensor is 

ImU* = tanh ( g ° + ( ^'^ )lm 77 Q/3 , (44) 

where n a p is the time ordered or causal polarization and is given by 

^jTr [T(G 2 ( P )J a G 4 (p + q)J )] . (45) 
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Above, /c M is the incoming neutrino four-momentum and is the four-momentum 
transfer. In writting the lepton tensor, we have neglected the electron mass term, 
since typical electron energies are of the order of a few hundred MeV. The Greens' 
functions Gi{p) (the index i labels particle species) describe the propagation of 
baryons at finite density and temperature. The current operator J M is 7^ for 
the vector current and 7^75 for the axial current. Effects of strong and elec- 
tromagnetic correlations may be calculated by utilizing the RPA polarization 
tensor 

jjRPA = n + JjRPA Dn ^ (46) 

where D denotes the interaction matrix, in Eq. (^2|) (see |47) for more details). 




Fig. 5. Neutrino mean free paths in matter with nucleons only (left panels). Right 
panels show ratios of mean free paths in matter without and with hyperons. Abscissa 
is baryon density n_g (no is the nuclear equilibrium density). Top panels show scatter- 
ing mean free paths common to all neutrino species. The bottom panels show results 
for electron neutrino mean free paths where absorption reactions are included. The 
neutrino content is labelled in the different panels 
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4.1 Neutrino Mean Free Paths 

The differential cross section (Eq. j42[) ) is needed in multi-energy group neutrino 
transport codes. However, more approximate neutrino transport algorithms (as 
in Section 2) require the total cross section as a function of the neutrino energy 
for the calculation of diffusion coefficients. The cross section per unit volume of 
matter (or equivalently the inverse mean free path) is obtained by integrating 
E 3 and n 3 in Eq. @. 

Under degenerate conditions even modest changes to the composition signif- 
icantly alter the neutrino scattering and absorption mean free paths. In Fig. || 
the neutrino scattering and absorption mean free paths are shown for models 
GM3np and GM3npH relevant to the deleptonization and cooling epochs. The 
top panels show the scattering mean free paths common to all neutrino species 
in neutrino free matter. The scattering mean free paths for thermal neutrinos 
(E u = 7rT) is shown in the left panel for various temperatures. To study the 
influence of hyperons, the ratio of the X np /X np H is shown in the right panels. 
The presence of hyperons significantly increase the scattering cross sections, by 
a factor ~ (2 — 3). Similar results for the absorption cross sections are shown in 
the lower panels for Yl = 0.4. Again we notice a significant enhancement (right 
panel) when hyperons appear, the factor here could be as large as 5. 



10.000 




Fig. 6. Left: Charged current inverse neutrino mean free paths versus temperature. 
Right: Comparison of scattering mean free paths in neutrino poor matter at fixed 
entropy for different EOSs in matter containing nucleons and also hyperons 

During the deleptonization stage, lepton number transport is sensitive to 
charged current reactions which dominate scattering reactions. At zero temper- 
ature, charged current reactions v + n <-> e + p depend sensitively on the proton 
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fraction Y p |58). Kinematic restrictions require Y p to be larger than 11 — 14% 
(direct Urea threshold). At early times, a finite neutrino chemical potential fa- 
vors a large Y p throughout the star, which enables these reactions to proceed 
without any hindrance. Toward the late stages, however, Y p decreases with de- 
creasing /i„ and charged current reactions may be naively expected to become 
inoperative. The threshold density for the charged current reaction when \i v = 
and T = depends sensitively on the density dependence of the nuclear symme- 
try energy. In field-theoretical models, in which the symmetry energy is largely 
given by contributions due to p-meson exchange, the critical density is typically 
iib = 2 ~ 3tiq. However, finite temperatures favor larger Y p 's and increase the 
average neutrino energy enabling the charged current reactions to proceed even 
below these densities. Fig. || shows that this is the case even at relatively low 
temperatures (T ~ 3 — 5) MeV for a baryon density ub = 0.15 fm~ 3 . The 
sharp rise with temperature, which occurs even for Y v = 0, clearly indicates that 
this reaction dominates the v e opacity even during the late deleptonization era. 
Thus, charged current reactions cannot be simply turned off when the neutrino 
chemical potential becomes small enough as was done in prior PNS simulations 

%■ 

The EOS and neutrino mean free paths are intimately related, which is best 
illustrated by comparing the results shown in Fig. || with those shown in Fig. ^. 
Composition and the baryon effective masses influence both the neutrino mean 
free paths and the matter's specific heat. Hyperons decrease the neutrino mean 
free paths at constant temperature Fig. |^. This trend is reversed at constant en- 
tropy due to the significantly lower temperatures favored in npH matter. Similar 
effects are apparent when we compare np models with different baryon effective 
masses. At a constant temperature, the larger effective mass in model GM3np 
favors larger cross sections, while at constant entropy this trend is again reversed 
due to the lower temperatures favored by the larger specific heat. 

The diffusion coefficients are calculated using Eq. ( |l8|) with the cross sections 
discussed above. The diffusion coefficients D2,Ds, and D4 are functions of Ub, 
T, and Y Vb . 



4.2 Inhomogeneous Phases: Effects of First Order Transitions 

The thermodynamics of the two situations, first order kaon condensation |l5| , ^| 
and the quark-hadron transition [^7| , has been previously considered. Reddy, 
Bertsch and Prakash |n| have studied the effects of inhomogeneous phases on 
^-matter interactions. Based on simple estimates of the surface tension between 
nuclear matter and the exotic phase, typical droplet sizes range from 5—15 fm 
|55) , and inter-droplet spacings range up to several times larger. The propagation 
of neutrinos whose wavelength is greater than the typical droplet size and less 
than the inter-droplet spacing, i.e., 2 MeV <J E v <J 40 MeV, will be greatly 
affected by the heterogeneity of the mixed phase, as a consequence of the coherent 
scattering of neutrinos from the matter in the droplet. 
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The Lagrangian that describes the neutral current coupling of neutrinos to 
the droplet is 

G 

where is the neutral current carried by the droplet and Gp — 1.166 x 
10~ 5 GeV -2 is the Fermi weak coupling constant. For non-relativistic droplets, 
= pw( x ) < ^ /i0 nas only a time like component. Here, pwi x ) is the ex- 
cess weak charge density in the droplet. The total weak charge enclosed in a 
droplet of radius r& is Nyy — J Q d d 3 x pw( x ) and the form factor is F(q) = 
(1/Nw) Jq" 1 d 3 x pw( x ) s'mqx/qx. The differential cross section for neutrinos 
scattering from an isolated droplet is then 

' /,T = ^G|<(l + cos0)F 2 ( <7 ). (48) 



dcos8 167T 

In the above equation, E v is the neutrino energy and 9 is the scattering angle. 
Since the droplets are massive, we consider only elastic scattering for which the 
magnitude of the momentum transfer is q — \/2E v (\ — cos 9). 

We must embed the droplets into the medium to evaluate the neutrino trans- 
port parameters. The droplet radius and the inter-droplet spacing are deter- 
mined by the interplay of surface and Coulomb energies. In the Wigner-Seitz 
approximation, the cell radius is Ryy = (3/4-ttNd) 1 ^ 3 where the droplet density 
is No- Except for one aspect, we will neglect coherent scattering from more 
than one droplet. If the droplets form a lattice, Bragg scattering will dominate 
and our description would not be valid. But for low density and a liquid phase, 
interference from multiple droplets affects scattering only at long wavelengths. 
If the ambient temperature is not small compared to the melting temperature, 
the droplet phase will be a liquid and interference effects arising from scatter- 
ing off different droplets are small for neutrino energies E u > (1/Rw)- However, 
multiple droplet scattering cannot be neglected for E v <J 1/R W . The effects of 
other droplets is to cancel scattering in the forward direction, because the inter- 
ference is destructive except at exactly zero degrees, where it produces a change 
in the index of refraction of the medium. These effects are usually incorporated 
by multiplying the differential cross section Eq. ( |48| ) by the static form factor of 
the medium. The static form factor, defined in terms of the radial distribution 
function of the droplets, g{r), is 

S{q) = 1 + N D J d 3 rexpiq.r (g(r) - 1) . (49) 

The droplet correlations, which determine g(r), arise due to the Coulomb force 
and is measured in terms of the dimensionless Coulomb number 
r = Z 2 e 2 1 (8nRwkT). Due to the long-range character of the Coulomb force, the 
role of screening and the finite droplet size, g(r) cannot be computed analytically. 
We use a simple ansatz for the radial distribution function g(r < Rw) — 
and g(r > Rw) = 1- For this choice, the structure factor is independent of r. 
Monte Carlo calculations [p9| of the liquid structure function of a simple one 
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component plasma indicate that our choice of S(q) is conservative for typical 
neutrino energies of interest. 

The simple ansatz for g(r) is equivalent to subtracting, from the weak charge 
density pw , a uniform density which has the same total weak charge Nw as the 
matter in the Wigner-Seitz cell. Thus, effects due to S(q) may be incorporated 
by replacing the form factor F(q) by 



F(q) -> F(q) = F(q) - 3 



sin qRw — (qRw) cos qRw 
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(50) 



The neutrino-droplet differential cross section per unit volume then follows: 
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Note that even for small droplet density Njj, the factor N w acts to enhance 
the droplet scattering. To quantify the importance of droplets as a source of 
opacity, we compare with the standard scenario in which matter is uniform and 
composed of neutrons. The dominant source of opacity is then due to scattering 
from thermal fluctuations and 
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where cy and ca are respectively the vector and axial coupling constants of the 
neutron, n n is the neutron number density, Ef n is the neutron Fermi energy and 
T is the matter temperature [ |60[ . 

The transport cross sections that are employed in studying the diffusive trans- 
port of neutrinos in the core of a neutron star arc differential cross sections 
weighted by the angular factor (1 — cos 9). The transport mean free path \(E V ) 
for a given neutrino energy E v is given by 
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(53) 



Models of first order phase transitions in dense matter provide the weak 
charge and form factors of the droplets and permit the evaluation of ^-droplet 
scattering contributions to the opacity of the mixed phase. For the models con- 
sidered, namely the first order kaon condensate and the quark-hadron phase 
transition, the neutrino mean free paths in the mixed phase are shown in the left 
and right panels of Fig. |^, respectively. The results are shown for the indicated 
values of the baryon density Ub and temperature T where the model predicts 
a mixed phase exists. The kaon droplets are characterized by radii rd ~ 7 fm 
and inter-droplet spacings Rw ~ 20 fm, and enclose a net weak vector charge 
Nw ~ 700. The quark droplets are characterized by r<2 ~ 5 fm and Rw ~ 11 
fm, and an enclosed weak charge Nw ~ 850. For comparison, the neutrino mean 
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Fig. 7. Neutrino mean free paths as a function of neutrino energy. Solid lines are for 
matter in a mixed phase containing kaons (left panel) and quarks (right panel), and 
dashed curves are for uniform matter. 



free path in uniform neutron matter at the same rib and T are also shown. It 
is apparent that there is a large coherent scattering-induced reduction in the 
mean free path for the typical energy E v ~ irT. At much lower energies, the 
inter-droplet correlations tend to screen the weak charge of the droplet, and at 
higher energies the coherence is attenuated by the droplet form factor. 

The large reduction in neutrino mean free path found here implies that the 
mixed phase will cool significantly slower than homogeneous matter. Conse- 
quently, the observable neutrino luminosity at late times might be affected as it 
is driven by the transport of energy from the deep interior. The reduced mean 
free path in the interior will tend to prolong the late time neutrino emission 
phase. 



4.3 Effects of Quark Superconductivity and Superfluidity 

Recent theoretical works |6l| , |62| suggest that quarks form Cooper pairs in medium, 
a natural consequence of attractive interactions destabilizing the Fermi surface. 
Although the idea of quark pairing in dense matter is not new |]6l| , p^| , it has 
recently seen renewed interest in the context of the phase diagram of QCD . 
Model calculations, mostly based on four-quark effective interactions, predict the 
restoration of spontaneously broken chiral symmetry through the onset of color 
superconductivity at low temperatures. They predict an energy gap of A ~ 100 
MeV for a typical quark chemical potential of \i q ~ 400 MeV. As in BCS theory, 
the gap will weaken for T > 0, and at some critical temperature T c there is 
a (second-order) transition to a "standard" quark-gluon plasma. During cool- 
ing from an initial temperature in excess of T c , the formation of a gap in the 
fcrmionic excitation spectrum in quark matter will influence various transport 
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properties of the system. Carter and Reddy have studied its influence on the 
transport of neutrinos p4| . 

The differential neutrino scattering cross section per unit volume in an in- 
finite and homogeneous system of relativistic fermions as calculated in linear 
response theory is given by Eq. (ff2|). The medium is characterized by the quark 
polarization tensor II a p. In the case of free quarks, each flavor contributes a 
term of the form 

n a p(q) = -iTr c J ^j4 Tr [ S o(p)r a S (p + q)T p } , (54) 

where So(?>) is the free quark propagator at finite chemical potential and temper- 
ature. The outer trace is over color and simplifies to a N c = 3 degeneracy. The 
inner trace is over spin, and the r a are the neutrino-quark vertex functions which 
determine the spin channel. Specifically, the vector polarization is computed by 
choosing [r a ,Tp) = (70,7/3). The axial and mixed vector- axial polarizations 
are similarly obtained from (r a ,rp) = (7075,7,375) and (fo,^) = (7^,7/375), 
respectively. 

The free quark propagators in Eq. (Q) are naturally modified in a supercon- 
ducting medium. As first pointed out by Bardccn, Cooper, and Schrieffer several 
decades ago, the quasi-particle dispersion relation is modified due to the presence 
of a gap in the excitation spectrum. In calculating these effects, we will consider 
the simplified case of QCD with two quark flavors which obey SU(2)^x SU(2)^ 
flavor symmetry, given that the light u and d quarks dominate low-energy phe- 
nomena. Furthermore we will assume that, through some unspecified effective 
interactions, quarks pair in a manner analogous to the BCS mechanism |65|] . The 
relevant consequences of this are the restoration of chiral symmetry (hence all 
quarks are approximately massless) and the existence of an energy gap at zero 
temperature, Aq, with approximate temperature dependence, 

A(T) = A 0] Jl-^ 2 . (55) 

The critical temperature T c ~ 0.57Z\o is likewise taken from BCS theory; this re- 
lation has been shown to hold for perturbative QCD |j(J and is thus a reasonable 
assumption for non-perturbative physics. 

Breaking the SU C (3) color group leads to complications not found in electro- 
dynamics. In QCD the superconducting gap is equivalent to a diquark conden- 
sate, which can at most involve two of the three fundamental quark colors. The 
condensate must therefore be colored. Since the scalar diquark (in the 3 color 
representation) appears to always be the most attractive channel, we consider 
the anomalous (or Gorkov) propagator J67| 

F{p)abfg = (qf a (p)Cj 5 q gb (-p)) 
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Here, a, b are color indices, /, g are flavor indices, e a b c is the usual anti-symmetric 
tensor and we have conventionally chosen 3 to be the condensate color. This 
propagator is also antisymmetric in flavor and spin, with C — —17072 being the 
charge conjugation operator. 

The color bias of the condensate forces a splitting of the normal quark propa- 
gator into colors transverse and parallel to the diquark. Quarks of color 3, parallel 
to the condensate in color space, will be unaffected and propagate freely, with 

S (p) b : f = ifi a S» f (^M, + (p„y - M70 ) . (57) 

\"o Po ^p J 

This is written in terms of the particle and anti-particle projection operators 
A + (p) and A~(p) respectively, where A (p) = (1 ± 707 • J>)/2. The excitation 
energies are simply E p = \p\ — fi for quarks and E p = \p\ + fj, for anti-quarks. 

On the other hand, transverse quark colors 1 and 2 participate in the diquark 
and thus their quasi-particle propagators are given as 

S(p)i 3 f = m + ^\ (p p y - M7o) • (58) 



The quasi-particle energy is £ p = y/ (\p\ ~ /j,) 2 + A 2 , and for the anti-particle 
4 = ^(\ P \ + f i) 2 + A 2 . 

The appearance of an anomalous propagator in the superconducting phase 
indicates that the polarization tensor gets contributions from both the normal 



quasi-particle propagators Eq. ( p8|) and anomalous propagator Eq. (56). Thus, 
to order G 2 F , Eq. ( |54| ) is replaced with the two contributions corresponding to 
the diagrams shown in Fig. §1 and written 



n a p{q) = -ij ?^i{ Tr [So{p)r a s Q {p + q)rp] 



+ 2Tr [S(p)r a S(p + q)r f3 ]+2Tr [F(p)r a F(p + q)r p ]} . (59) 

The remaining trace is over spin, as the color trace has been performed. Fig. ||(a) 
corresponds to the first two terms, which have been decomposed into one term 
with ungapped propagators Eq. ( |57| ) and the other with gapped quasi-particle 
propagators Eq. ([58]). Fig. ||(b) represents the third, anomalous term. 

For neutrino scattering we must consider vector, axial, and mixed vector-axial 
channels, all summed over flavors. The full polarization, to be used in evaluating 
Eq. (p2[), may be written 

n aP = ]T [{clfnl p + {c{fn^ - 2C f v c{n^] . (60) 



/ 



The coupling constants for up quarks are Cy = i — | sin 9w and C\ = ^ , and 
for down quarks, C v = — \ + | sin 2 8 W and C\ = — |, where sin 2 6 W ~ 0.23 is 
the Weinberg angle. 
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Fig. 9. Left panel: Neutrino mean free path as a function of neutrino energy E v . Right 
panel: Neutrino mean free paths for E v = nT as a function of A/T. These results are 
virtually independent of temperature for T <J 50 MeV. 



The differential cross section, Eq. ([42]) and the total cross section are ob- 
tained by integrating over all neutrino energy transfers and angles. Results for 
the neutrino mean free path, A = V/cr, are shown in Fig. |^ as a function of 
incoming neutrino energy E v (for ambient conditions of \x q = 400 MeV and 
T = 30 MeV). They show the same energy dependence found previously for free 
relativistic and degenerate fermionic matter [Q; A oc 1/E% for E v ^> T and 
A oc 1/E V for E v -C T. The results indicate that this energy dependence is not 
modified by the presence of a gap when A ~ T. Thus, the primary effect of 
the superconducting phase is a much larger mean free path. This is consistent 
with the suppression found in the vector-longitudinal response function, which 
dominates the polarization sum Eq. (|60|), at go < q. 

We now consider the cooling of a macroscopic sphere of quark matter, a 
toy approximation for the core of a neutron star with a mixed quark phase, as 
it becomes superconducting. As in the preceding calculation, we consider the 
relatively simple case of two massless flavors with identical chemical potentials 
and disregard the quarks parallel in color to the condensate; i.e. we consider a 
background comprised exclusively of quasi-quarks. 
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The cooling of a spherical system of quark matter from T ~ T c ~ 50 MeV 
is driven by neutrino diffusion, for the neutrino mean free path is much smaller 
than the dimensions of system of astrophysical size, and yet several orders of 
magnitude larger than the mean free path of the quarks. The diffusion equation 
for energy transport by neutrinos in a spherical geometry is 



dT 
~dt 



1 8L V 
r 2 dr 



(61) 



where Cy is the specific heat per unit volume of quark matter, T is the tem- 
perature, and r is the radius. The neutrino energy luminosity for each neutrino 
type, L v , depends on the neutrino mean free path and the spatial gradients in 
temperature and is approximated by an integral over neutrino energy E v 



L„ = 6 J dE u -^E 3 v r 2 X(E v ) 



df{E v ) 



dr 



(62) 



We assume that neutrino interactions are dominated by the neutral current 
scattering which is common to all neutrino types, accounting for the factor 6 in 
Eq. ©. 

The solution to the diffusion equation will depend on the initial temperature 
gradients. However, being primarily interested in a qualitative description of 
cooling through a second-order phase transition to superconducting matter, the 
temporal behavior can be characterized by a time scale t c which is proportional 
to the inverse cooling rate. This characteristic time is 



t c (T) = C V (T) 



R 2 



c(A(T)) ' 



(63) 



and is a strong function of the temperature since it depends on both the specific 
heat Cy and the energy-averaged neutrino mean free path, (A(T)). The latter 
quantity is here approximated by \(E V — ttT) since the neutrinos are in thermal 
equilibrium. Using BCS theory, as described in the previous section, (A(T)) 
depends on the gap A as shown in the right panel of Fig. ^. The results indicate 
that for small A/T, A is not strongly modified, but as A/T increases so too does 
A, non-linearly at first and then exponentially for A/T> 5. Also, in the BCS 
theory, the temperature dependences of Cy (dashed curve) and A (dot-dashed 
curve) are shown in Fig. [10[ Finally, the ratio t^(T)/t c (T), a measure of the 
extent to which the cooling rate is changed by a gap, is shown by the solid line 
in Fig. [l(| We note that the diffusion approximation is only valid when A <C R 
and will thus fail for very low temperatures, when A <J R. 

These results are readily interpreted. The cooling rate around T c is influenced 
mainly by the peak in the specific heat associated with the second order phase 
transition, since the neutrino mean free path is not strongly affected when T > 
A. Subsequently, as the matter cools, both Cy and A -1 decrease in a non-linear 
fashion for T ~ A. Upon further cooling, when T <C A, both Cy and A -1 
decrease exponentially. Both effects accelerate the cooling process. 
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Fig. 10. The extent to which different physical quantities are affected by a supercon- 
ducting transition. Ratios of the cooling time scale (solid curve), the inverse mean 
free path (short-dashed curve) and the matter specific heat (dot-dashed curve) in the 
superconducting phase to that in the normal phase is shown as a function of the mat- 
ter temperature. The ratio of the gap to its zero temperature value A is also shown 
(long-dashed curve). 

We conclude that if it were possible to measure the neutrino luminosity 
from the hypothetical object described here, a second order superconducting 
transition might be identified by the temporal characteristics of the late time 
supernova neutrino signal from a PNS. Specifically, there might be a brief interval 
during which the cooling would slow when the core temperature falls below T c , 
signified by a period of reduced neutrino detection. However, this effect might be 
obscured by v— opaque matter outside the star's core. If the neutrino opacity of 
these outer regions of the star is large, it is likely that any sharp temporal feature 
associated with neutrino transport in the core will be diluted as the neutrino 
diffuse through the outer regions. Nevertheless, the main finding, which is that 
phase transitions in the core can have a discernible impact on the transport of 
neutrinos and suggests that the late time supernova neutrino signal is a promising 
probe of the high density and low temperature region of the QCD phase diagram. 

5 Results of PNS Simulations 

Neutrino signals from PNSs depends on many stellar properties, including the 
mass; initial entropy, lepton fraction and density profiles; and neutrino opaci- 
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ties. Pons et al. |T(J carried out a detailed study of the dependence of neutrino 
emission on PNS characteristics. They verified the generic results of Burrows & 
Lattimer jlj that both neutrino luminosities and average energies increased with 
increasing mass (see Fig. |IT| ). In addition, they found that variations in initial 
entropy and lepton fraction profiles in the outer regions of the PNS caused only 
transient (lasting a few tenths of a second) variations in neutrino luminosities 
and energies. Variations in the central lepton fraction and entropy were found 
to produce modest changes in neutrino luminosities that persisted to late times. 
The central values of lepton fraction and entropy are established during core 
collapse, and will depend upon the initial properties of the star as well as the 
EOS and neutrino transport during the collapse. 




Fig. 11. Left: The evolution of the neutrino average energy and total neutrino lumi- 
nosity is compared for several assumed PNS baryon masses and EOSs. The EOSs in the 
left panels contain only baryons and leptons while those in the right panels also contain 
hyperons. Top right panel: Evolution of the central baryon number density for different 
baryonic mass stars containing hyperons (model GM3npH) which are metastable. Bot- 
tom right panel: Time required by stars shown in the top panel to reach the unstable 
configuration 



5.1 Baseline Results 

Properties of the dense matter EOS that affect PNS evolution include the com- 
pressibility symmetry energy, specific heat, and composition. Pons et al. |l0| 
employed a field theoretical EOS (J], with which the results due to some differ- 
ences in stellar size (due to variations in nuclear interactions) and composition 
were studied. Some results are summarized in Fig. O. Overall, both the average 
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energies and luminosities of stars containing hyperons were larger compared to 
those without. In addition, for stars without hyperons, those stars with smaller 
radii had higher average emitted neutrino energy, although the predicted lu- 
minosities for early times (t < 10 s) were insensitive to radii. This result only 
holds if the opacities are calculated consistently with each EOS [E6p7j]; other- 
wise rather larger variations in evolutions would have been found ]42| , |48| . The 
same held true for models which allowed for the presence of hyperons, except 
when the initial proto-neutron star mass was significantly larger than the maxi- 
mum mass for cold, catalyzed matter. Another new result was that the average 
emitted neutrino energy of all flavors increased during the first 2-5 seconds of 
evolution, and then decreased nearly linearly with time. For times larger than 
about 10 seconds, and prior to the occurrence of neutrino transparency, the 
neutrino luminosities decayed exponentially with a time constant that was sen- 
sitive to the high-density properties of matter. Significant variations in neutrino 
emission occurred beyond 10 seconds: it was found that neutrino luminosities 
were larger during this time for stars with smaller radii and with the inclusion 
of hyperons in the matter. Finally, significant regions of the stars appeared to 
become convectively unstable during the evolution , as several works have found 

The right panel of Fig. [ll| shows the time development of the central baryon 
density (top panel) and also the time to the collapse instability as a function 
of baryon mass (bottom panel). The larger the mass, the shorter the time to 
instability, since the PNS does not have to evolve in lepton number as much. 
Above 2.005 M Q , the metatstability disappears because the GM3npH initial 
model with the lepton and entropy profiles we chose is already unstable. Below 
about 1.73 Mq, there is no metastability, since this is the maximum mass of the 
cold, catalyzed npH star for GM3. The signature of neutrino emission from a 
metastable PNS should be identifiable and it is discussed in Section 5.4. 



5.2 Influence of Many-Body Correlations 

The main effect of the larger mean free paths produced by RPA corrections |l7J 
is that the inner core deleptonizes more quickly. In turn, the maxima in central 
temperature and entropy are reached on shorter timescales. In addition, the 
faster increase in thermal pressure in the core slows the compression associated 
with the deleptonization stage, although after 10 s the net compressions of all 
models converge. 

The relatively large, early, changes in the central thermodynamic variables 
do not, however, translate into similarly large effects on observables such as the 
total neutrino luminosity and the average radiated neutrino energy, relative to 
the baseline simulation. The luminosities for the different models are shown as 
a function of time in Fig. [l^. The left panel shows the early time development 
in detail. The exploratory models agree with the results reported in ]4^^9|| . 
However, the magnitude of the effects when full RPA corrections are applied is 
somewhat reduced compared to the exploratory models. It is especially impor- 
tant that at and below nuclear density, the corrections due to correlations are 
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Fig. 12. Left: The upper panel shows the total emitted neutrino luminosity for the 
PNS evolutions described in Reddy et al. p7| |. The lower panel shows the ratio of 
the luminosities obtained in the three models which contain corrections to the baseline 
(Hartree approximation) model. Right: Emitted neutrino luminosity for long-term PNS 
evolutions 



relatively small. Since information from the inner core is transmitted only by 
the neutrinos, the time scale to propagate any high density effect to the neu- 
trinosphere is the neutrino diffusion time scale. Since the neutrinosphere is at a 
density approximately 1/100 of nuclear density, and large correlation corrections 
occur only above 1/3 nuclear density where nuclei disappear, we find that corre- 
lation corrections calculated here have an effect at the neutrinosphere only after 
1.5 s. Moreover, the RPA suppresion we have calculated is considerably smaller 
than those reported in [4^^9| , reaching a maximum of about 30% after 5 s, com- 
pared to a luminosity increase of 50% after only 2 s. However, the corrections 
are still very important during the longer-term cooling stage (see Fig. |l^), and 
result in a more rapid onset of neutrino transparency compared to the Hartree 
results. 



5.3 Signals in Detectors 

In Fig. [l^ the lifetimes versus M b for stars containing hyperons (npH) and 
npK stars are compared In both cases, the larger the mass, the shorter 
the lifetime. For kaon-rich PNSs, however, the collapse is delayed until the final 
stage of the Kelvin-Hclmholtz epoch, while this is not necessarily the case for 
hyperon-rich stars. 

In Fig. |l4| the evolution of the total neutrino energy luminosity is shown 
for different models. Notice that the drop in the luminosity for the stable star 
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(solid line), associated with the end of the Kelvin- Hclmholtz epoch, occurs at 
approximately the same time as for the metastable stars with somewhat higher 
masses. In all cases, the total luminosity at the end of the simulations is below 
10 51 erg/s. The two upper shaded bands correspond to SN 1987A detection limits 
with KII and 1MB, and the lower bands correspond to detection limits in SNO 
and SuperK for a future galactic supernova at a distance of 8.5 kpc. The times 
when these limits intersect the model luminosities indicate the approximate times 
at which the count rate drops below the background rate (dN / 'dt) bg = 0.2 Hz. 

The poor statistics in the case of SN 1987A precluded a precise estimate 
of the PNS mass. Nevertheless, had a collapse to a black hole occurred in this 
case, it must have happened after the detection of neutrinos ended. Thus the 
SN 1987A signal is compatible with a late kaonization-induced collapse, as well 
as a collapse due to hyperonization or to the formation of a quark core. More 
information would be extracted from the detection of a galactic SN with the new 
generation of neutrino detectors. 
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Fig. 13. Lifetimes of metastable stars as a function of the stellar baryon mass. Solid 
lines show results for PNSs containing kaon-condensates and dashed lines show the 
results of Pons, et al. [[LCfl for PNSs containing hyperons 



In SNO, about 400 counts are expected for electron antineutrinos from a 
supernova located at 8.5 kpc. The statistics would therefore be improved signif- 
icantly compared to the observations of SN 1987A. A sufficiently massive PNS 
with a kaon condensate becomes metastable, and the neutrino signal terminates, 
before the signal decreases below the assumed background. In SuperK, however, 
up to 6000 events are expected for the same conditions (because of the larger 
fiducial mass) and the effects of metastability due to condensate formation in 
lower mass stars would be observable. 
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Fig. 14. The evolution of the total neutrino luminosity for stars of various baryon 
masses. Shaded bands illustrate the limiting luminosities corresponding to a count 
rate of 0.2 Hz in all detectors, assuming a supernova distance of 50 kpc for 1MB and 
Kamioka, and 8.5 kpc for SNO and SuperK. The width of the shaded regions represents 
uncertainties in the average neutrino energy from the use of a diffusion scheme for 
neutrino transport 

5.4 What Can We Learn Prom Neutrino Detections? 

The calculations of Pons, et al. [^3) show that the variations in the neutrino 
light curves caused by the appearance of a kaon condensate in a stable star are 
small, and are apparently insensitive to large variations in the opacities assumed 
for them. Relative to a star containing only nucleons, the expected signal differs 
by an amount that is easily masked by an assumed PNS mass difference of 
0.01 — 0.02 M Q . This is in spite of the fact that, in some cases, a first order 
phase transition appears at the star's center. The manifestations of this phase 
transition are minimized because of the long neutrino diffusion times in the star's 
core and the Gibbs' character of the transition. Both act in tandem to prevent 
either a "core-quake" or a secondary neutrino burst from occurring during the 
Kelvin-Helmholtz epoch. 

Observable signals of kaon condensation occur only in the case of metastable 
stars that collapse to a black hole. In this case, the neutrino signal for a star 
closer than about 10 kpc is expected to suddenly stop at a level well above that 
of the background in a sufficiently massive detector with a low energy threshold 
such as SuperK. This is in contrast to the signal for a normal star of similar mass 
for which the signal continues to fall until it is obscured by the background. The 
lifetime of kaon-condensed metastable stars has a relatively small range, of order 
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50-70 s for the models studied here, which is in sharp contrast to the case of 
hyperon-rich metastable stars for which a significantly larger variation in the 
lifetime (a few to over 100 s) was found. This feature of kaon condensation 
suggests that stars that destabilize rapidly cannot do so because of kaons. 

Pons, et al. |7^| determined the minimum lifetime for metastable stars with 
kaons to be about 40 s by examining the most favorable case for kaon conden- 
sation, which is obtained by maximizing the magnitude of the optical potential. 
The maximum optical potential is limited by the binary pulsar mass constraint, 
which limits the star's maximum gravitational mass to a minimum value of 1.44 
Mq. Therefore, should the neutrino signal from a future supernova abruptly ter- 
minate sooner than 40 s after the birth of the PNS, it would be more consistent 
with a hyperon- or quark-induced instability than one due to kaon condensation. 

It is important to note that the collapse to a black hole in the case of kaon 
condensation is delayed until the final stages of the Kelvin-Helmholtz epoch, 
due to the large neutrino diffusion time in the inner core. Consequently, to 
distinguish between stable and metastable kaon-rich stars through observations 
of a cessation of a neutrino signal from a galactic supernovae is only possible 
using sufficiently massive neutrino detectors with low energy thresholds and low 
backgrounds, such as the current SNO and SuperK, and future planned detectors 
including the UNO. 



5.5 Expectations From Quark Matter 

Strangeness appearing in the form of a mixed phase of strange quark matter also 
leads to metastability. Although quark matter is also suppressed by trapped neu- 
trinos H?],^], the transition to quark matter can occur at lower densities than 
the most optimistic kaon case, and the dependence of the threshold density on 
Yl is less steep than that for kaons. Thus, it is an expectation that metasta- 
bility due to the appearance of quarks, as for the case of hyperons, might be 
able to occur relatively quickly. Steiner, et al. |22| have demonstrated that the 
temperature along adiabats in the quark-hadron mixed phase is much smaller 
than what is found for the kaon condensate-hadron mixed phase. Calculations 
of PNS evolution with a mixed phase of quark matter, including the possible 



effects of quark matter superfluidity 64 arc currently in progress and will be 
reported separately. 



6 Long-Term Cooling: The Next Million Years 

Following the transparency of the neutron star to neutrinos, the only observa- 
tional link with these objects is through photon emissions, either as a pulsar or 
through thermal emissions or both. Thermal emissions of course are controlled 
by the temperature evolution of the star, and this depends sensitively upon 
its internal composition. The tabulation of temperatures and ages for a set of 
neutron stars would go a long way to deciding among several possibilities. 
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6.1 Thermal Evolution 

The cooling of a young (age < 10 5 yr) neutron star is mainly governed by 
v— emission processes and the specific heat j74|. Due to the extremely high ther- 
mal conductivity of electrons, a neutron star becomes nearly isothermal within 
a time t w w 1 — 100 years after its birth, depending upon the thickness of the 
crust Q. After this time its thermal evolution is controlled by energy balance: 

f-ftf-^-UI, (64) 

where E t h is the total thermal energy and Gy is the specific heat. L 7 and L v 
are the total luminosities of photons from the hot surface and vs from the in- 
terior, respectively. Possible internal heating sources, due, for example, to the 
decay of the magnetic field or friction from differential rotation, are included in 
H. Our cooling simulations were performed by solving the heat transport and 
hydrostatic equations including general relativistic effects (see [[r||). The sur- 
face's effective temperature T e is much lower than the internal temperature T 
because of a strong temperature gradient in the envelope. Above the envelope 
lies the atmosphere where the emerging flux is shaped into the observed spec- 
trum from which T e can be deduced. As a rule of thumb T e /10 6 K w ^/T/10 8 K, 
but modifications due to magnetic fields and chemical composition may occur. 



6.2 Rapid vs. Slow Cooling 

The simplest possible v emitting processes are the direct Urea processes fa + 
£ — > fa + ve , fa — > /i + t + ~t, where fa and fa are either baryons or quarks 
and £ is either an electron or a muon. These processes can occur whenever 
momentum conservation is satisfied among fa, fa and £ (within minutes of birth, 
the v chemical potential vanishes). If the unsuppressed direct Urea process for 
any component occurs, a neutron star will rapidly cool because of enhanced 
emission: the star's interior temperature T will drop below 10 9 K in minutes 
and reach 10 7 K in about a hundred years. T e will hence drop to less than 
300,000 K after the crustal diffusion time t w [0j7|Q. This is the so-called 
rapid cooling paradigm. If no direct Urea processes are allowed, or they are all 
suppressed, cooling instead proceeds through the significantly less rapid modified 
Urea process in which an additional fermion enables momentum conservation. 
This situation could occur if no hyperons are present, or the nuclear symmetry 
energy has a weak density dependence |5£|[77]]. The v emisssion rates for the 
nuclcon, hyperon, and quark Urea and modified Urea processes can be found 




6.3 Superfluid and Superconducting Gaps 

Pairing is unavoidable in a degenerate Fermi liquid if there is an attractive 
interaction in any channel. The resulting superfluidity, and in the case of charged 
particles, superconductivity, in neutron star interiors has a major effect on the 
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star's thermal evolution through suppressions of neutrino (y) emission processes 
and specific heats 0,^5) ■ Neutron (ra) , proton (p) and /1-hyperon superfluidity in 
the 1 Sq channel and n superfluidity in the 3 P2 channel have been shown to occur 
with gaps of a few MeV or less |7^,Q. However, the density ranges in which gaps 
occur remain uncertain. At large baryon densities for which perturbative QCD 
applies, pairing gaps for like quarks have been estimated to be a few MeV p3| |. 
However, the pairing gaps of unlike quarks (ud, us, and ds) have been suggested 
to be several tens to hundreds of MeV through non-perturbative studies [j62| 
kindling interest in quark superfluidity and superconductivity |8l||82|| and their 
effects on neutron stars [}S3[|84|| . 

The effect of the pairing gaps on the emissivities and specific heats for massive 
baryons are investigated in Q and are here generalized to the case of quarks. 
The principal effects are severe suppressions of both the emissivity and specific 
heat when T << A, where A is the pairing gap. In a system in which several su- 
perfluid species exist the most relevant gap for these suppressions is the smallest 
one. The specific heat suppression is never complete, however, because leptons 
remain unpaired. Below the critical temperature T c , pairs may recombine, re- 
sulting in the emission of vv pairs with a rate that exceeds the modified Urea 
rate below 10 10 K p6| ]; these processes are included in our calculations. 

The baryon and quark pairing gaps we adopt are shown in Fig. |l5| Note 
that gaps are functions of Fermi momenta i denoting the species) which 

translates into a density dependence. For pp(n,p) <; 200 — 300 MeV/c, nucleons 
pair in the 1 So state, but these momenta correspond to densities too low for 
enhanced v emission involving nucleons to occur. At higher p_f's, baryons pair 
in higher partial waves. The n 3 P2 gap has been calculated for the Argonne 
Vis, CD-Bonn and Nijmegen I & II interactions 79 . This gap is crucial since it 
extends to large p_p(n) and can reasonably be expected to occur at the centers 
of neutron stars. For > 350 MeV/c, gaps are largely uncertain because 

of both experimental and theoretical uncertainties |79| . The curves [a] , [b] and 
[c] in Fig. |l5| reflect the range of uncertainty. The p 3 P2 gap is too small to be 
of interest. Gaps for the 1 So pairing of A, taken from Q and shown as dotted 
curves, are highly relevant since As participate in direct Urea emission as soon 
as they appear j7?J. Experimental information beyond the 1 So channel for A 
is not available. As for S— hyperons remain largely unexplored. The quark (q) 
gaps are taken to be Gaussians centered at pf{<1) — 400 MeV/c with widths of 
200 MeV/c and heights of 100 MeV [model D], 10 MeV [C], 1 MeV [B] and 0.1 
MeV [A] , respectively. The reason for considering quark gaps much smaller than 
suggested in [ |62"[|65| is associated with the multicomponent nature of charge- 
neutral, beta-equilibrated, neutron star matter as will become clear shortly. 



6.4 Effects of Composition 

We consider four generic compositions: charge-neutral, beta equilibrated mat- 
ter containing nucleons only (np), nucleons with quark matter (npQ), nucleons 
and hyperons (npH), and nucleons, hyperons and quarks (npHQ). In the cases 
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Fig. 15. Pairing gaps adopted for neutron 1 So and 3 P2, proton 1 So, A 1 So, and quarks. 
The n 3 P2 gaps are anisotropic; plotted values are angle- averaged. The A gaps corre- 
spond, in order of increasing A, to background densities ns = 0.48, 0.64 and 0.8 fm~ 3 , 
respectively. The s-wave quark gaps are schematic; see text for details 



involving quarks, a mixed phase of baryons and quarks is constructed by satis- 
fying Gibbs' phase rules for mechanical, chemical and thermal equilibrium p^ |. 
The phase of pure quark matter exists only for very large baryon densities, and 
rarely occurs in our neutron star models. Baryonic matter is calculated using a 
field-theoretic model at the mean field level ; quark matter is calculated us- 
ing either a bag- like model or the Nambu- Jona-Lasinio quark model | 
equation of state (EOS) is little affected by the pairing phenomenon, since the 
energy density gained is negligible compared to the ground state energy densities 
without pairing. 

Additional particles, such as quarks or hyperons, have the effect of softening 
the EOS and increasing the central densities of stars relative to the np case. 
For the npQ model studied, a mixed phase appears at the density ns — 0.48 
fm -3 . Although the volume fraction of quarks is initially zero, the quarks them- 
selves have a significant pf(q) when the phase appears. The pfS of the three 
quark flavors become the same at extremely high density, but for the densities 
of interest they are different due to the presence of negatively charged leptons. 
In particular, pf(s) is much smaller than pf(u) and pp(d) due to the larger 
s-quark mass. Use of the Nambu-Jona-Lasinio model, in which quarks acquire 
density-dependent masses resembling those of constituent quarks, exaggerates 
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the reduction of pf(s). This has dramatic consequences since the pairing phe- 
nomenon operates at its maximum strength when the Fermi momenta are exactly 
equal; even small asymmetries cause pairing gaps to be severely reduced (8^j8^| . 
In addition, one may also expect p-wave superfluidity, to date unexplored, which 
may yield gaps smaller than that for the s-wave. We therefore investigate pairing 
gaps that are much smaller than those reported for the case of s-wave superflu- 
idity and equal quark pf's. 

The introduction of hyperons does not change these generic trends. In the 
case npH, the appearance of hyperons changes the lepton and nucleon pf's 
similarly to the appearance of quarks although with less magnitude. While the 
appearance of quarks is delayed by the existence of hyperons, at high densities 
the pf's of nucleons and quarks remain similar to those of the npQ case. The 
existence of either hyperons or quarks, however, does allow the possibility of 
additional direct Urea processes involving themselves as well as those involving 
nucleons by decreasing pj?(n) — Pf(p)- For the npQ and npHQ models studied, 
the maximum masses are = 1.5M© , the central baryon densities are = 1.35 
fm~ 3 , and the volume fractions of quarks at the center are = 0.4. 



6.5 Examples of Results 

Cooling simulations of stars without hyperons and with hyperons are compared, 
in Figs. [l6] and 17, respectively, to available observations of thermal emissions 



from pulsars. Sources for the observational data can be found in [ p9| . However, at 
the present time, the inferred temperatures must be considered as upper limits 
because the total flux is contaminated, and in some cases dominated, by the pul- 
sar's magnetospheric emission and/or the emission of a surrounding synchrotron 
nebula. Furthermore, the neutron star surface may be reheated by magneto- 
spheric high energy photons and particles; late-time accretion for non-pulsing 
neutron stars is also possible. Other uncertainties arise in the temperature esti- 
mates due to the unknown chemical composition and magnetic field strength in 
the surface layers, and in the age, which is based upon the observed spin-down 
time. In these figures, the bolder the data symbol the better the data. 



np and npQ Matter 

The np case is considered in Fig. |l6[ in which solid lines indicate the temperature 
evolution of a 1.4 M Q star for quarkless matter: case [z] is for no nucleon pairing 
at all, and cases [a] , [b] and Jcl correspond to increasing values for the neutron 
3 P2 gap, according to Fig. [15]. The field-theoretical model employed for the 
nucleon interactions allows the direct nucleon Urea process, which dominates the 
cooling. The unimpeded direct Urea process carries the temperature to values 
well below the inferred data. Pairing suppresses the cooling for T < T c , where T 
is the interior temperature, so T e increases with increasing A. If the direct Urea 
process is not allowed, the range of predicted temperatures is relatively narrow 
due to the low cmissivity of the modified Urea process. We show an example of 




Fig. 16. Cooling of 1.4MQ stars with np matter (continuous curves) and npQ matter 
(dashed and dotted curves). The curves labelled as [a], [b], and [c] correspond to n 3 P2 
gaps as in Fig. [z] corresponds to zero n gap. Models labelled [A], [B], [C] and [D] 
correspond to quark gaps as in Fig. |l^; [Z] corresponds to zero quark gap 



such cooling (curve [Mc]) using the n 3 P2 gap [c] for a 1.4M Q with an EOS |50| 
for which the direct Urea cooling is not allowed. 

The other curves in the figure illustrate the effects of quarks upon the cooling. 
The dotted curves [Z] are for vanishingly small quark gaps; the dashed curves 



([A], [B], [C] and [D]) are for quark gaps as proposed in Fig. 15. For nonexistent 
([z]) or small ([a]) nucleon gaps, the quark Urea process is irrelevant and the de- 
pendence on the existence or the size of the quark gaps is very small. However, for 
large nucleon gaps ([b] and [c]), the quark direct Urea process quickly dominates 
the cooling as the nucleon direct Urea process is quenched. It is clear that for 
quark gaps of order 1 MeV or greater ([B], [C] or [D]) the effect of quarks is again 
very small. There is at most a slight increase in the stars temperatures at ages 
between 10 1 to 10 5 ~ 6 years due to the reduction of pF(n) and the consequent 
slightly larger gap (Fig. |l^). Even if the quark gap is quite small ([A]), quarks 
have an effect only if the nucleon gap is very large ([b] or [c]), i.e., significantly 
larger than the quark gap: the nucleon direct Urea process is suppressed at high 
temperatures and the quark direct Urea process has a chance to contribute to 
the cooling. We find that the effects of changing the stellar mass M are similar 
to those produced by varying the baryon gap, so that only combinations of M 
and A might be constrained by observation. 
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Fig. 17. Cooling of stars with npH (continuous lines) and npHQ matter (dotted lines) 
for various stellar masses (in Mq ). n 3 P2 gaps are from case [c] while quark gaps, when 
present, are from model [C] of Fig. [l^ 



npH and npHQ Matter 



The thermal evolution of stars containing hyperons has been discussed in 19 




but we obtain qualitatively different results here. Hyperons open new direct Urea 
channels: A — > p + e + V e and A + e — * S~ + v e if I7 _ 's are present, with their 
inverse processes. Previous results showed that the cooling is naturally controlled 
by the smaller of the A or n gap. However, this is significantly modified when 
the A gap is more accurately treated. At the A appearance threshold, the gap 
must vanish since pf(A) is vanishingly small. We find that a very thin layer, 
only a few meters thick, of unpaired or weakly paired A's is sufficient to control 
the cooling. This effect was overlooked in previous works perhaps because they 
lacked adequate zonal resolution. 

In Fig. [ll] we compare the evolution of stars of different masses made of 
either npH or npHQ matter. We find that all stars, except the most massive 
npH ones, follow two distinctive trajectories depending on whether their central 
density is below or above the threshold for A appearance (= 0.54 fm -3 in our 
model EOS, the threshold star mass being 1.28Mq ). In the case of npH matter, 
stars with M > I.50M Q are dense enough so that the A 1 So gap vanishes and 
hence undergo fast cooling, while stars made of npHQ matter do not attain such 
high densities. The temperatures of npH stars with masses between 1.3 and 1.5 
M Q are below the ones obtained in the models of Fig. |l6| with the same n 3 P2 
gap [b], which confirms that the cooling is dominated by the very thin layer 
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of unpaired /l's (the slopes of these cooling curves are typical of direct Urea 
processes). Only if the n 3 P2 gap ^0.3 MeV do the cooling curves fall below 
what is shown in Fig. |l7]. Notice, moreover, that in the mass range 1.3 - 1.48 
Mq the cooling curves are practically indistinguishible from those with unpaired 
quark matter shown in Fig. [l^. In these models with npH or npHQ matter, 
there is almost no freedom to "fine-tune" the size of the gaps to attain a given 
T e : stars with yl's will all follow the same cooling trajectory, determined by the 
existence of a layer of unpaired or weakly paired A's, as long as the n 3 P2 gap 
is not smaller. It is, in some sense, the same result as in the case of np and npQ 
matter: the smallest gap controls the cooling and now the control depends on 
how fast the A 1 So gap increases with increasing pp(A). 

6.6 Implications 

Our results indicate that observations could constrain combinations of the smaller 
of the neutron and A— hyperon pairing gaps and the star's mass. Deducing the 
sizes of quark gaps from observations of neutron star cooling will be extremely 
difficult. Large quark gaps render quark matter practically invisible, while van- 
ishing quark gaps lead to cooling behaviors which are nearly indistinguishable 
from those of stars containing nucleons or hyperons. Moreover, it also appears 
that cooling observations by themselves will not provide definitive evidence for 
the existence of quark matter itself. 

7 The Structure of Catalyzed Stars 

In this section, we explore from a theoretical perspective, how the structure of 
neutron star depends upon the assumed EOS. This study is crucial if new ob- 
servations of masses and radii are to lead to effective constraints of the EOS 
of dense matter. Two general classes of stars can be identified: normal stars in 
which the density vansihes at the stellar surface, and self-bound stars in which 
the density at the surface is finite. Normal stars originate from nuclear force 
models which can be conveniently grouped into three broad categories: non- 
relativistic potential models, relativistic field theoretical models, and relativistic 
Dirac-Brucckncr-Hartree-Fock models. In each of these approaches, the presence 
of additional softening components such as hyperons, Bose condensates or quark 
matter, can be incorporated. Details of these approaches have been further con- 
sidered in Lattimer et al. J93| and Prakash et al. Q. A representative sample, 
and some general attributes, including references and typical compositions, of 
equations of state employed here are summarized in Table |l|. 

For normal matter, the EOS is that of interacting nucleons above a transition 
density of 1/3 to 1/2 n s . Below this density, the ground state of matter consists 
of heavy nuclei in equilibrium with a neutron-rich, low-density gas of nucleons. 
However, for most of the purposes of this paper, the pressure in the region 
n < 0.1 fm -3 is not relevant as it does not significantly affect the mass-radius 
relation or other global aspects of the star's structure. Nevertheless, the value 
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of the transition density, and the pressure there, arc important ingredients for 
the determination of the size of the superfluid crust of a neutron star that is 
believed to be involved in the phenomenon of pulsar glitches (Link, Epstein, & 
Lattimer @). 

Four equations of state are taken from Akmal & Pandharipande Jll[]. These 
are: API (the AV18 potential), AP2 (the AV18 potential plus Svb relativistic 
boost corrections), AP3 (the AV18 potential plus a three-body UIX potential 
), and AP4 (the AV18 potential plus the UIX potential plus the Svb boost). 
Three equations of state from Muller & Serot |j(J, labelled MS1-3, correspond 
to different choices of the parameters £ and £ which determine the strength of 
the nonlinear vector and isovector interactions at high densities. The numer- 
ical values used are £ = C = 0;£ = 1.5, £ = 0.06; and £ = 1.5, £ = 0.02, 
respectively. Six EOSs come from the phcnomcnological non-relativistic poten- 
tial model of Prakash, Ainsworth & Lattimer p5|, labelled PAL1-6, which have 
different choices of the symmetry energy parameter at the saturation density, its 
density dependence, and the bulk nuclear matter incompressibility parameter 
K s . The incompressibilities of PAL1-5 were chosen to be K s = 180 or 240 MeV, 
but PAL6 has K s = 120 MeV. Three interactions denoted GM1-3 come from 
the field-theoretical model of Glendenning & Moszkowski |Q . Two interactions 
come from the field-theoretical model of Glendenning & Schaffner-Bielich JDS] : 
GL78 with U K (p ) = -140 MeV and TM1 with U K = -185 MeV. The la- 
bels denoting the other EOSs in Table [l] are identical to those in the original 
references. 

The rationale for exploring a wide variety of EOSs, even some that are rela- 
tively outdated or in which systematic improvements are performed, is two-fold. 
First, it provides contrasts among widely different theoretical paradigms. Second, 
it illuminates general relationships that exist between the pressure-density rela- 
tion and the macroscopic properties of the star such as the radius. For example, 
AP4 represents the most complete study to date of Akmal & Pandharipande 
[ [Tl"| , in which many-body and special relativistic corrections are progressively 
incorporated into prior models, AP1-3. AP1-3 are included here because they 
represent different pressure-energy density-baryon density relations, and serve 
to reinforce correlations between neutron star structure and microscopic physics 
observed using alternative theoretical paradigms. Similarly, several different pa- 
rameter sets for other EOSs are chosen. 

In all cases, except for PS (Pandharipande & Smith p7|), the pressure is 
evaluated assuming zero temperature and beta equilibrium without trapped neu- 
trinos. PS only contains neutrons among the baryons, there being no charged 
components. We chose to include this EOS, despite the fact that it has been su- 
perceded by more sophisticated calculations by Pandharipande and coworkers, 
because it represents an extreme case producing large radii neutron stars. 

The pressure-density relations for some of the selected EOSs are shown in 
Fig. |l8| which displays three significant features to note for normal EOSs. First, 
there is a fairly wide range of predicted pressures for beta-stable matter in the 
density domain n s /2 < n < 2n s . For the EOSs displayed, the range of pressures 
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Table 1. Approach refers to the underlying theoretical technique. Composition refers 
to strongly interacting components (n=neutron, p=proton, H=hyperon, K=kaon, 
Q=quark); all models include leptonic contributions. 



Symbol 



Reference 




Approach 


Composition 


Variational 


np 


Potential 


mr 


Variational 


np 


Variational 


np 


Field Theoretical 


np 


Dirac-Brueckner HF 


np 


Dirac-Brueckner HF 


np 


Schematic Potential 


np 


Field Theoretical 


npH 


Field Theoretical 


npK 


Field Theoretical 


npHQ 


Quark Matter 


Q(u, d, s) 



FP 
PS 

WFF(l-3) 

AP(l-4) 

MS(l-3) 

MPA(l-2) 

ENG 

PAL (1-6) 

GM(l-3) 

GS(l-2) 

PCL(l-2) 

SQM(l-3) 



Friedman & Pandharipande 
Pandharipande & Smith 
Wiringa, Fiks & Fabrocine 
Akmal & Pandharipande [til 
Miiller & Serot 

Miither, Prakash & Ainsworth 



Engvik et al. |10C ] 
Prakash, Ainsworth & Lattimer |9fJ 
Glendenning & Moszkowski jl3j 
Glendenning & Schaffner-Bielich J^|] 



Prakash, Cooke & Lattimer 
Prakash, Cooke & Lattimer 
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covers about a factor of five, but this survey is by no means exhaustive. That 
such a wide range in pressures is found is somewhat surprising, given that each of 
the EOSs provides acceptable fits to experimentally-determined nuclear matter 
properties. Clearly, the extrapolation of the pressure from symmetric matter to 
nearly pure neutron matter is poorly constrained. Second, the slopes of the pres- 
sure curves are rather similar. A polytropic index of n ~ 1, where P — Kn 1+1 / n , 
is implied. Third, in the density domain below 2n s , the pressure-density relations 
seem to fall into two groups. The higer pressure group is primarily composed of 
relativistic field-theoretical models, while the lower pressure group is primarily 
composed of non-relativistic potential models. It is significant that relativistic 
field-theoretical models generally have symmetry energies that increase propor- 
tionately to the density while potential models have much less steeply rising 
symmetry energies. 

A few of the plotted normal EOSs have considerable softening at high densi- 
ties, especially PAL6, GS1, GS2, GM3, PS and PCL2. PAL6 has an abnormally 
small value of incompressibility (K s — 120 MeV). GS1 and GS2 have phase tran- 
sitions to matter containing a kaon condensate, GM3 has a large population of 
hyperons appearing at high density, PS has a phase transition to a neutral pion 
condensate and a neutron solid, and PCL2 has a phase transition to a mixed 
phase containing strange quark matter. These examples are representative of the 
kinds of softening that could occur at high densities. 

The best-known example o f sel f-bound stars results from Wit ten's [ L0 1 1 con- 
jecture (also see Fahri & Jaffc 102 , Haensel, Zdunik & Schaeffcr [ 103 1 , Alcock & 
Olinto [ 104 1 , and Prakash et al. |105|| ) that strange quark matter is the ultimate 
ground state of matter. In this paper, the self-bound EOSs are represented by 
strange-quark matter models SQM1-3, using perturbative QCD and an MIT- 
type bag model, with parameter values given in Table g. The existence of an 
energy ceiling equal to the baryon mass, 939 MeV, for zero pressure matter 
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Fig. 18. Left: The pressure-density relation for a selected set of EOSs contained in 
Table [l]. The pressure is in units of MeV fm~ 3 and the density is in units of baryons 
per cubic fermi. The nuclear saturation density is approximately 0.16 fm -3 . Right: The 
pressures inferred from the empirical correlation Eq. (^), for three hypothetical radius 
values (10, 12.5 and 15 km) overlaid on the pressure-density relations shown on the 
left. The light shaded region takes into account only the uncertainty associated with 
C(n, M); the dark shaded region also includes a hypothetical uncertainty of 0.5 km in 
the radius measurement. The neutron star mass was assumed to be 1.4 Mq 

requires that the bag constant B < 94.92 MeV fm~ 3 . This limiting value is 
chosen, together with zero strange quark mass and no interactions (a c = 0), for 
the model SQM1. The other two models chosen, SQM2 and SQM3, have bag 
constants adjusted so that their energy ceilings are also 939 MeV. 

Table 2. Parameters for self-bound strange quark s tars . Numerical values employed 
in the MIT bag model as described in Fahri & Jaffe |l02| . 



Model 


B (MeV fin -3 ) 


m s (MeV) 




SQM1 


94.92 








SQM2 


64.21 


150 


0.3 


SQM3 


57.39 


50 


0.6 
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7.1 Neutron Star Radii 

Fig. |l9| displays the resulting mass-radius relations for catalyzed matter. Rhoades 
] demonstrated that the assumption of causality beyond a fidu- 
cial density pf sets an upper limit to the maximum mass of a neutron star: 
4.2-^/ Ps/pf M Q . Lattimer et al. J93| have shown that the causality constraint 
also sets a lower limit to the radius: R JsS 1.52i? s , where R s = 2GM/c 2 , which 
is shown in Fig. 19. For a 1.4 M© star, this is about 4.5 km. The most reliable 
estimates of neutron star radii in the near future will likely stem from observtins 
of thermal emission from their surfaces. Such estimates yield the so-called "radi- 
ation radius" i?oo = Rf^l — R s , a quantity resulting from redshifting the star's 
luminosity and temperature. A given value of R^ imiplies that R < R^ and 
M < 0.13(i?oo/km) M Q . Contours of Roo are also displayed. With the exception 
of model GS1, the EOSs used to generate Fig. [l^ result in maximum masses 
greater than 1.442 M Q , the limit obtained from PSR 1913+16. From a theo- 
retical perspective, it appears that values of R^ in the range of 12-20 km are 
possible for normal neutron stars whose masses are greater than 1 M . 

One observes that normal neutron stars have minimum masses of about 0.1 
M that are primarily determined by the EOS below n s . At the minimum mass, 
the radii are generally in excess of 100 km. Self-bound stars have no minimum 
mass and the maximum mass self-bound stars have nearly the largest radii pos- 
sible for a given EOS. If the strange quark mass m s = and interactions are 
neglected (a c = 0), the maximum mass is related to the bag constant B in the 
MIT-type bag model by M max = 2.033 (56 MeV far 3 / 5 ) 1/2 M o- Prakash et al. 
| fl05| and Lattimer et al. |33) showed that the addition of a finite strange quark 
mass and/or interactions produces larger maximum masses. The constraint that 
M max > 1.44 Mq is thus automatically satisfied by the condition that the energy 
ceiling is 939 MeV, and non-zero values of m s and a c yield larger radii for every 
mass. The locus of maximum masses is given simply by R = 1.85R S (Lattimer 
et al. |9^]) as shown in the right-hand panel of Fig. |l9[ Strange quark stars wit h 
electrostatically supported normal-matter crusts (Glendenning & Weber |107| |) 
have larger radii than those with bare surfaces. Coupled with the additional 
constraint M > 1M Q from protoneutron star models, MIT-model strange quark 
stars cannot have R < 8.5 km or R^ < 10.5 km. These values are comparable 
to the smallest possible radii for a Bose (pion or kaon) condensate EOS. 

One striking feature of Fig. is that in the mass range from 1-1.5 M Q or 
more the radius has relatively little dependence upon the stellar mass. The major 
exceptions illustrated are the model GS1, in which a mixed phase containing a 
kaon condensate appears at a relatively low density, and the model PAL6, which 
has an extremely small nuclear incompressibility (120 MeV). Both of these have 
considerable softening and a large increase in central density for M > 1 M . 

While it is generally assumed that a stiff EOS implies both a large maximum 
mass and a large radius, many counter examples exist. For example, GM3, MSI 
and PS have relatively small maximum masses but have large radii compared to 
most other EOSs with larger maximum masses. Also, not all EOSs with extreme 
softening have small radii for M > 1 M Q (e.g., GS2, PS). Nonetheless, for stars 
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Fig. 19. Mass-radius curves for several EOSs listed in Table |l|. The left panel is for 
stars containing nucleons and, in some cases, hyperons. The right panel is for stars 
containing more exotic components, such as mixed phases with kaon condensates or 
strange quark matter, or pure strange quark matter stars. In both panels, the lower 
limit causality places on R is shown as a dashed line, a constraint derived from glitches 
in the Vela pulsar is shown as the solid line labelled AI/I — 0.014, and contours of 
constant Roc = R/ yjl — 2GM/Rc 2 are shown as dotted curves. In the right panel, 
the theoretical trajectory of maximum masses and radii for pure strange quark matter 
stars is marked by the dot-dash curve labelled R — 1.85R S 

with masses greater than 1 M Q , only models with a large degree of softening 
(including strange quark matter configurations) can have Roo < 12 km. 

To understand the relative insensitivity of the radius to the mass for normal 
neutron stars, it is relevant that a Newtonian polytrope with n = 1 has the 
property that the stellar radius is independent of both the mass and central 
density. An n = 1 polytrope also has the property that the radius is proportional 
to the square root of the constant K in the polytropic pressure law P = Kp 1+1 / n . 
This suggests that there might be a quantitative relation between the radius and 
the pressure that does not depend upon the EOS at the highest densities, which 
determines the overall softness or stiffness (and hence, the maximum mass). 

In fact, this conjecture may be verified. Fig. shows a remarkable empirical 
correlation between the radii of 1 and 1.4 Mq normal stars and the matter's 
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Fig. 20. Empirical relation between pressure, in units of MeV fm -3 , and 7?, in km, for 
EOSs listed in Table [[]. The upper panel shows results for 1 M Q (gravitational mass) 
stars; the lower panel is for 1.4 Mq stars. The different symbols show values of i?P _1//4 
evaluated at three fiducial densities 



pressure evaluated at fiducial densities of 1, 1.5 and 2 n s . Numerically, the cor- 
relation has the form of a power law between the radius Rm, defined as the 
radius at a particular mass M, and the total pressure P(n) evaluated at a given 
density: 

Rm C(n, M) [P(7i)] a23 -°- 26 . (65) 

C(n, M) is a number that depends on the density n at which the pressure was 
evaluated and the stellar mass M. An exponent of 1/4 was chosen for display 
in Fig. ^o], but the correlation holds for a small range of exponents about this 
value. Using an exponent of 1/4, and ignoring points associated with EOSs with 
phase transitions in the density ranges of interest, we find values for C(n, M), in 
units of km fm 3 / 4 MeV -1 / 4 , which are listed in Table § The error bars are taken 
from the standard deviations. The correlation is seen to be somewhat tighter for 
the baryon density n — 1.5n s and 2n s cases. Despite the relative insensitivity 
of radius to mass for a particular EOS in this mass range, the nominal radius 
Rm has a variation ~ 5 km. The largest deviations from this correlation occur 
for EOSs with extreme softening or for configurations very near their maximum 
mass. This correlation is valid only for cold, catalyzed neutron stars, i.e., not 
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for protoneutron stars which have finite entropies and might contain trapped 
neutrinos. 

Table 3. The quantity C(n,M) of Eq. (||) The quantity C(n,M), in units of km 
fm 3/4 MeV _:L/4 , which relates the pressure (evaluated at density n) to the radius of 
neutron stars of mass M. The errors are standard deviations 



n 


1 Mq 


1.4 M 


n s 


9.53 ±0.32 


9.30 ± 0.60 


1.5n s 


7.14 ±0.15 


7.00 ±0.31 


2n s 


5.82 ±0.21 


5.72 ±0.25 



If a measurement of P near n s can be deduced in this way, an important 
clue about the symmetry properties of matter will be revealed. The energy per 
particle and pressure of cold, beta stable nucleonic matter is 

E(n, x) ~ E(n, 1/2) + S v (n)(l - 2xf , 

P{n,x) ~n 2 [E'(n,l/2) + S' v (n)(l-2x) 2 }, (66) 

where E(n, 1/2) is the energy per particle of symmetric matter and S v (n) is the 
bulk symmetry energy (which is density dependent). Primes denote derivatives 
with respect to density. If only one term in this expansion is important, as noted 
by Prakash, Ainsworth & Lattimer jojj, then 

S v (n) ~ \ d * E ^ X) ~ E ^ °) - E{n > V 2 ) ■ ( 6? ) 

At n a , the symmetry energy can be estimated from nuclear mass systematics and 
has the value S v = S v (n s ) ~ 27—36 MeV. Attempts to further restrict this range 
from consideration of fission barriers and the energies of giant resonances provide 
constraints between S v and S v (n) primarily by providing correlations between 
S v and S s , the surface symmetry parameter. Lattimer & Prakash 1 108 detail 
how S s is basically a volume integral of the quantity 1 — S v /S v (n) through the 
nucleus. However, both the magnitude of S v and its density dependence S v (n) 
remain uncertain. Part of the bulk symmetry energy is due to the kinetic energy 
for noninteracting matter, which for degenerate nucleonic matter is proportional 
to n 2 / 3 , but the remainder of the symmetry energy, due to interactions, is also 
expected to contribute significantly to the overall density dependence. 

Leptonic contributions must to be added to Eq. ( |66|) to obtain the total en- 
ergy and pressure; the electron energy per baryon is (3/4)hcx(^TT 2 nx) 1 ^ 3 . Matter 
in neutron stars is in beta equilibrium, i.e., fj, e = fi n — fi p = —dE/dx, which 
permits the evaluation of the equilibrium proton fraction. The pressure at the 
saturation density becomes 



P s =n s (l- 2x a ) [n s S' v (l - 2x s ) + S v x s ] , 



(68) 
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where S' v = S' v (n s ) and the equilibrium proton fraction at n s is 

x s ~ {3TT 2 n s )- 1 (AS v /hc) 3 ~ 0.04, (69) 

for S„ = 30 MeV. Due to the small value of x s , we find that P s ~ n 2 s S' v . 

Were we to evaluate the pressure at a larger density, contributions fea- 
turing other nuclear parameters, including the nuclear incompressibility K s = 
9(dP/dn)\ ns and the skewness K' s = —27nl(d 3 E/dn 3 )\ ns , also contribute. How- 
ever, the K s and K' s terms largely cancel, up to 2n s , so the symmetry term 
dominates. 

At present, experimental guidance concerning the density dependence of the 
symmetry energy is limited and mostly based upon the division of the nu- 
clear symmetry energy between volume and surface contributions. Upcoming 
experiments involving heavy-ion collisions which might sample densities up to 
~ (3 — 4)n s , will be limited to analyzing properties of the nearly symmetric 
nuclear matter EOS through a study of matter, momentum, and energy flow of 



nucleons. However, the parity-violating experiment [ 109 1 to accurately determine 
the thickness of the neutron skin in 208 Pb at Jefferson Lab will provide impor- 
tant constraints. The neutron skin thickness is directly proportional to S s / S v . In 
addition, studies of heavy nuclei far off the neutron drip lines using radioactive 
ion beams might also provide useful constraints. 



7.2 Moments of Inertia 



Besides the stellar radius, other global attributes of neutron stars are potentially 
observable, including the moment of inertia and the binding energy. These quan- 
tities depend primarily upon the ratio M/R as opposed to details of the EOS 
(Lattimer & Prakash ||). 

The moment of inertia, for a star uniformly rotating with angular velocity 
Q, is 

7=(87r/3)/ r 4 (p + P/c 2 )e<- x -^/ 2 (Lu/f2)dr. (70) 

(71) 



The metric function uj{t) is a solution of the equation 

d[r A e- {x+v)/2 J]/dr + A^ude^+^l 2 / dr 
with the surface boundary condition 



R 

ujr ; = Q uj'n = fl 

3 



1 - 



2GI 



(72) 



The second equality in the above follows from the definition of / and the TOV 
equation. Writing j = exp[— (y + A)/2], the TOV equation becomes 



f = -4:TTGr(P/c 2 + p)je x /c 2 



Then, one has 



2^ 
3G 



r 3 dj 



,3 p4, ,/ 



(73) 
(74) 
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Useful approximations which are valid for three analytic, exact, solutions to 
GR, the incompressible fluid (Inc), the Tolman VII (Tolman JH0[; VII) solution, 



and BuchdahPs [111] solution (Buch), are 



I Inc /MR 2 ~ 2(1 -0.87/3 - 0.3/3 2 )- 1 / 5 , (75) 
iBuch/MR 2 ~ (2/3 - 4/^ 2 )(l - 1.81/3 + 0.47/3 2 )" 1 , (76) 
Itvii/MR 2 ~ 2(1- 1.1/3 -0.6/3 2 )"V7. (77) 

Fig. |2l] indicates that the T VII approximation is a rather good approximation to 
most EOSs without extreme softening at high densities, for M/R > 0.1 M Q /km. 
The EOSs with softening fall below this trajectory. 
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Fig. 21. The moment of inertia /, in units of MR 2 , for several EOSs listed in Table |l|. 
The curves l abell ed "Inc", "T VII", "Buch" a nd " RP" are for an incompressible fluid, 
the Tolman |llO[ VII solution, the Buchdahl [111] solution, and an approximation of 
Ravenhall & Pethick [112], respectively. The inset shows details of I /M R 2 for M /R — > 



Another interesting result from Fig. ^T] concerns the moments of inertia of 
strange quark matter stars. Such stars are relatively closely approximated by 
incompressible fluids, this behavior becoming exact in the limit of /3 — > 0. This 
could have been anticipated from the M oc i? 3 behavior of the M—R trajectories 
for small (3 strange quark matter stars as observed in Fig. n9l 
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7.3 Crustal Fraction of the Moment of Inertia 

A new observational constraint involving / concerns pulsar glitches. Occasionally, 
the spin rate of a pulsar will suddenly increase (by about a part in 10 6 ) without 
warning after years of almost perfectly predictable behavior. However, Link, 



Epstein & Lattimer 94 argue that these glitches are not completely random: 
the Vela pulsar experiences a sudden spinup about every three years, before 
returning to its normal rate of slowing. Also, the size of a glitch seems correlated 
with the interval since the previous glitch, indicating that they represent self- 
regulating instabilities for which the star prepares over a waiting time. The 
angular momentum requirements of glitches in Vela imply that > 1.4% of the 
star's moment of inertia drives these events. 

Glitches are thought to represent angular momentum transfer between the 
crust and another component of the star. In this picture, as a neutron star's 
crust spins down under magnetic torque, differential rotation develops between 
the stellar crust and this component. The more rapidly rotating component 
then acts as an angular momentum reservoir which occasionally exerts a spin- 
up torque on the crust as a consequence of an instability. A popular notion at 
present is that the freely spinning component is a superfluid flowing through a 
rigid matrix in the thin crust, the region in which dripped neutrons coexist with 
nuclei, of the star. As the solid portion is slowed by electromagnetic forces, the 
liquid continues to rotate at a constant speed, just as superfluid He continues 
to spin long after its container has stopped. This superfluid is usually assumed 
to locate in the star's crust, which thus must contain at least 1.4% of the star's 
moment of inertia. 

The high-density boundary of the crust is naturally set by the phase boundary 
between nuclei and uniform matter, where the pressure is Pt and the density 
n t . The low-density boundary is the neutron drip density, or for all practical 
purposes, simply the star's surface since the amount of mass between the neutron 
drip point and the surface is negligible. One can utilize Eq. ( f70| ) to determine 
the moment of inertia of the crust alone with the assumptions that P j c 2 < < p, 
m(r) ~ M, and toj ~ u>r and P oc p 4 / 3 in the crust (Lattimer & Prakash H: 



Al ^ 287rP t i? 3 (1 - 1.67/? - 0.6/3 2 ) 
T ~ 3Mc 2 ]) 

In general, the EOS parameter P t , in the units of MeV fm -3 , varies over the 
range 0.25 < Pt < 0.65 for realistic EOSs. The determination of this parameter 
requires a calculation of the structure of matter containing nuclei just below 
nuclear matter density that is consistent with the assumed nuclear matter EOS. 
Unfortunately, few such calculations have been performed. Like the fiducial pres- 
sure at and above nuclear density which appears in Eq. (|6^), Pt should depend 
sensitively upon the behavior of the symmetry energy near nuclear density. 

Link, Epstein & Lattimer established a lower limit to the radius of the 
Vela pulsar by using Eq. ( |78| ) with P t at its maximum value and the glitch 



2P t (l + 5/3- 14/3 



2\- ' 



n t mbC 2 (3 2 
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Fig. 11. Mass-radius curves for selected EOSs from Table [j], comparing theoretical 
contours of AI/I — 0.014 from appr oximations developed in this paper, labelled "LP", 
and from Ravenhall & Pethick [112], labelled "RP", to numerical results (solid dots). 



Two values of Pt, the transition pressure demarking the crust's inner boundary, which 
bracket estimates in the literature, are employed. The region to the left of the P t = 0.65 
MeV fm -3 curve is forbidden if Vela glitches are due to angular momentum transfers 
between the crust and core, as discussed in Link, Epstein & Lattimer [Q. For compar- 
ison, the region excluded by causality alone lies to the left of the da shed curve labelled 
"causality" as determined by Lattimer et al. 193] and Glendenning 113] 



constraint AI/I > 0.014: 

R > 3.9 + 3.5A//M - O.O8(M/M ) 2 km . (79) 

As shown in Fig. ^2], this constraint is somewhat more stringent than one based 
upon causality. Better estimates of the maximum value of Pt should make this 
constraint more stringent. 



7.4 Binding Energies 

The binding energy formally represents the energy gained by assembling N 
baryons. If the baryon mass is mj, the binding energy is simply BE = Nrnt — AI 
in mass units. However, the quantity mb has various interpretations in the lit- 
erature. Some authors take it to be 939 MeV/c 2 , the same as the neutron or 
proton mass. Others take it to be about 930 MeV/c 2 , corresponding to the mass 
of C 12 /12 or Fe 56 /56. The latter choice would be more appropriate if BE was to 
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represent the energy released in by the collapse of a white-dwarf-like iron core in 
a supernova explosion. The difference in these definitions, 10 MeV per baryon, 
corresponds to a shift of fO/939 ~ 0.01 in the value of BE/M . This energy, 
BE, can be deduced from neutrinos detected from a supernova event; indeed, it 
might be the most precisely determined aspect of the neutrino signal. 



Lattimer & Yahil 114] suggested that the binding energy could be approxi- 
mated as 

BE w 1.5 • 10 51 (M/M Q ) 2 ergs = O.O84(M/M ) 2 M Q . (80) 

Prakash et al. Q also concluded that such a formula was a reasonable approxi- 
mation, based upon a comparison of selected non-relativistic potential and field- 
theoretical models, good to about ±20 %. 

However, Lattimer & Prakash Q proposed a more accurate representation 
of the binding energy: 

BE/M ~ 0.6/3/(1 - 0.5/3) , (81) 

which incorporates some radius dependence. Thus, the observation of supernova 
neutrinos, and the estimate of the total radiated neutrino energy, will yield more 
accurate information about M/R than about M alone. 

7.5 Outlook for Radius Determinations 

Any measurement of a radius will have some intrinsic uncertainty. In addition, 
the empirical relation we have determined between the pressure and radius has 
a small uncertainty. It is useful to display how accurately the equation of state 
might be established from an eventual radius measurement. This can be done 
by inverting Eq. (|6^), which yields 

P(n)~[R M /C(n,M)} A . (82) 

The inferred ranges of pressures, as a function of density and for three possible 
values of Ria, are shown in the right panel of Fig. [l^. It is assumed that the 
mass is 1.4 Mq, but the results are relatively insensitive to the actual mass. Note 
from Table ||that the differences between C for 1 and 1.4 Mq are typically less 
than the errors in C itself. The light shaded areas show the pressures including 
only errors associated with C. The dark shaded areas show the pressures when 
a hypothetical observational error of 0.5 km is also taken into account. These 
results suggest that a useful restriction to the equation of state is possible if the 
radius of a neutron star can be measured to an accuracy better than about 1 
km. 

The reason useful constraints might be obtained from just a single measure- 
ment of a neutron star radius, rather than requiring a series of simultaneous 



mass-radii measurements as Lindblom |115 proposed, stems from the fact that 



we have been able to establish the empirical correlation, Eq. fl65|). In turn, it 
appears that this correlation exists because most equations of state have slopes 
d In P/d In n ~ 2 near n s . 
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8 Tasks and Prospects 

There are several topics that will merit attention from theoretical and experi- 
mental perspectives. Among those dealing with i/-matter interactions are: 

Dynamic Structure Functions from Microscopic Calculations 

Neutrino properties of astrophysical interest depend crucially on the nature of 
the excitation spectrum of the nuclear medium to spin and spin-isospin probes. 
The excited states are very different depending upon whether or not interac- 
tions in the medium conserve spin and spin-isospin. The importance of tensor 
correlations in the medium is thus clear, since they break the conservation laws. 



Friman & Maxwell [116] hrst emphasized the importance of tensor correlations 
in the process 

v e + n + n — ► e~ + p + n , (83) 

and noted that their neglect underestimates the rate of v e absorption by as much 
as an order of magnitude. In their study, they used a hard core description of 
the short range correlations and a one pion exchange model for the medium and 



long range ones. Sawyer & Soni [117 and Haensel & Jerzak |118{ , who used 
additional correlations based on a Reid soft core potential, confirmed that large 
reductions were possible in degenerate matter for non-degenerate neutrinos. 

Since these earlier works, many-body calculations have vastly improved (e.g. 
[pi]) and have been well-tested against data on light nuclei and nuclear matter. 
Much better tensor correlations are now available, so that we may better pin 
down the rate of absorption due to the above process. Detailed calculations 
to include arbitrary matter and neutrino degeneracies encountered in in many 
astrophysical applications are necessary. 

Axial Charge Renormalization 



In dense matter, the axial charge of the baryons is renormalized [119 120 ,121 1 , 
which alters the neutrino-baryon couplings from their vacuum values. Since the 
axial contribution to the scattering and absorption reactions is typically three 
times larger than the vector contributions, small changes in the axial vector 
coupling constants significantly affect the cross sections. The calculation of this 
renormalization requires a theoretical approach which treats the pion and chiral 
symmetry breaking explicitly. So far, this has been done in isospin symmetric 



nuclear matter [122], but not for neutron matter or for beta-equilibrated neu- 
tron star matter. Substantial reductions may be expected in the ^-matter cross 
sections from this in- medium effect. 



Multi-Pair Excitations 

Neutrinos can also excite many-particle states in an interacting system, inverse 
bremsstrahlung being an example of a two-particle excitation [ 123 1 . These excita- 
tions provide an efficient means of transferring energy between the neutrinos and 
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baryons which are potentially significant in low-density matter. However, multi- 
group neutrino transport will be needed to include this effect. In addition, such 
calculations require source terms for neutrino processes such as bremsstrahlung 
and neutrino pair production. The latter process has been accurately treated in 




7r and K Dispersion Relations Through //-Nucleus Reactions 

The experimental program that would do the most to illuminate theoretical 
issues permeating neutrino interactions in dense matter would be studies of 
neutrino reactions on heavy nuclei, the only direct way of probing the matrix 
elements of the axial current in nuclear matter. Pioneering suggestions in this 



regard have been put forth by Sawyer & Soni [125 126[ , Ericson [127 , and Sawyer 
[ p.28| . The basic idea is to detect positively charged leptons or e + ) produced 
in inclusive experiments 

v + X^n + ( ore+) ( orK")+X (84) 

which is kinematically made possible when the in-medium 7r~ or K~ dispersion 
relation finds support in space-like regions. The sharp peaks at forward angles 
in the differential cross section versus lepton momentum survive the 100-200 
MeV width in the incoming GeV or so neutrinos from accelerator experiments. 
Calculations of the background from quasi-elastic reactions indicate that the 
signal would be easily detectable. 



i/-Matter Interactions at Sub-Nuclear Density 

Analgous to the effects of inhomogeneities for j/-matter interactions discussed 
in the case of a first-order kaon condensate or quark matter transition is the 
case of coherent scattering of neutrinos from closely-spaced nuclei at sub-nuclear 
densities. The sizes and separations of nuclei are similar to those of the droplets 
discussed for the kaon and quark situations, so the range of neutrino energies 
most affected will be similar. These will be important in reshaping the v spectrum 
from PNSs, and are of potential importance in the supernova mechanism itself 
due to the large energy dependence of {/-matter cross sections behind the shock. 

In addition to these, several topics of interest from an astrophysical perspec- 
tive include: 



Improvements in PNS Simulations 

These include a) an adequate treatment of convection coupled with neutrino 
transport appear to be necessary based upon large regions that are potentially 
convectively unstable; b) the consideration of other softening components in 
dense matter that might produce effects dissimilar to those found when consid- 
ering hyperons and kaons, i.e., quarks; c) improved transport calculations with 
many energy groups, especially in the transparent regime; d) a self-consistent 



Evolution of a Neutron Star 



57 



treatment of accretion, which is known to significantly contribute to the early v 
emission. The latter two items necessitate the coupling of a multi-group trans- 
port scheme with a hydrodynamical code of the type generally used for supernova 
simulations. 



Determination of the Neutron Skin of Neutron-rich Nuclei 

The Jefferson Lab experiment ]109{ (PREX) is anticipated to yield accurate 
measurements of the neutron-skin thickness of 208 Pb. This quantity, from a the- 
oretical viewpoint, is the volume integral of the inverse of the symmetry energy 
throughout a nucleus, and represents how the nuclear symmetry energy is split 
between volume and surface contributions. Not coincidentally, the density de- 
pendence of the symmetry energy is also implicated in the predicted neutron 
star radius M. 



Determination of the Radius of a Neutron Star 

The best prospect for measuring a neutron star's radius may be the nearby 



object RX J185635-3754. Parallax information [129] indicates its distance to be 
about 60 pc. In addition, it may be possible to identify spectral lines with the 
Chandra and XMM X-ray facilities that would not only yield the gravitational 
redshift, but would identify the atmospheric composition. Not only would this 
additional information reduce the uncertainty in the deduced value of i?oo, but, 
both the mass and radius for this object might thereby be estimated. It is also 
possible that an estimate of the surface gravity of the star can be found from 
further comparisons of observations with atmospheric modelling, and this would 
provide a further check on the mass and radius. 
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